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ABSTRACT 

We have used the MHW2 filter to obtain estimates of the flux densities at the WMAP frequencies 
of a complete sample of 2491 sources, mostly brighter than 500 mJy at 5 GHz, distributed over 
the whole sky excluding a strip around the Galactic equator < 5°). After having detected 933 
k> | sources at the > 3er level in the MHW2 filtered maps - our New Extragalactic WMAP Point Source 

(NEWPS) 3ct Catalogue - we are left with 381 sources at > 5<t in at least one WMAP channel, 369 
of which constitute our NEWPS 5(T catalogue. It is remarkable to note that 98 (i.e. 26%) sources 
detected at > 5er are 'new" , they are not present in the WMAP catalogue. Source fluxes have been 
corrected for the Eddington bias. Our flux density estimates before such correction are generally 
in good agreement with the WMAP ones at 23 GHz. At higher frequencies WMAP fluxes tend to 
be slightly higher than ours, probably because WMAP estimates neglect the deviations of the point 
spread function from a Gaussian shape. On the whole, above the estimated completeness limit of 1.1 
Jy at 23 GHz we detected 43 sources missed by the blind method adopted by the WMAP team. On 
the other hand, our low-frequency selection threshold left out 25 WMAP sources, only 12 of which, 
however, are > 5cr detections and only 3 have S23GHZ > 1-1 Jy Thus, our approach proved to be 
competitive with, and complementary to the WMAP one. 
Subject headings: filters: wavelets; point sources: catalogues, identifications 
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1. INTRODUCTION 

As a by-product of its temperature and polarization 
maps of the Cosmic Microwave Background (CMB), the 
Wilkinson Microwave Anisotropy Probe (WMAP) mis- 
sion has produced the first all-sky surveys of extragalac- 
tic sources at 23, 33, 41, 61 and 94 GHz (Bennett et al. 
2003; Hinshaw et al. 2006), a still unexplored frequency 
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region where many interesting astrophysical phenomena 
are expected to show up (see e.g. De Zotti et al. 2005). 

^From the analysis of the first three years survey data 
the WMAP team has obtained a catalogue of 323 ex- 
tragalactic point sources (EPS; Hinshaw et al. 2006), 
substantially enlarging the first-year one that included 
208 EPS detected above a flux limit of ~ 0.8-1 Jy 
(|Bennett et al.ll2003Q. 

As discussed by [Hinshaw et al.l (|2006l ). the detection 
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process - quite similar to the one adopted for producing 
the first-year catalogue - can be summarized as follows. 

1/2 

First, the weighted map, N a ^ s T, where N b s is the num- 
ber of observations per pixel, is filtered in the harmonic 
space by the global matched filter, 

bf , „ 
(1) 

^2^fcmb _j_ ^fnoisc ' V ' 

where bg is the transfer function of the WM AP beam 
response (|Page et al.ll2fM Uarosik et al.|[2006h . C| mb is 
the angular power spectrum of the CMB and C" OISC rep- 
resents the noise power. Then, all the peaks above the 5a 
threshold in the filtered maps are assumed to be source 
detections and are fitted in real space, i.e. in the unfil- 
tered maps, to a Gaussian profile plus a planar baseline. 
The Gaussian amplitude is finally converted to a flux 
density using the conversion factors given in Table 5 of 
iPage et al.l(|2003D . The error on the flux density is given 
by the statistical error of the fit. 

If a source has been detected above 5a, a being the 
noise rms defined globally, in a frequency channel, the 
flux densities in the other channels are given if they are 
above 2a and the source width falls within a factor of two 
of the true beam width. In other words, not all sources 
listed at any give n frequency in the th ree-year WMAP 
catalogue of EPS (jHinshaw et al.ll2006l Table 9) are 5a 
detections at that frequency. Finally, to identify the 
detected sources with 5 GHz counterparts, they cross- 
correl ated them with the GB6, PMN and iKuhr et al.1 
(1981) catalogues. 

The fact that almost all sources detected by WMAP 
were previously catalogued at lower frequencies suggests 
that a fuller exploitation of the WMAP data can be 
achieved complementing the blind search already carried 
out by the WMAP team with a search for WMAP coun- 
terparts to sources detected at lower frequencies. The 
latter approach exploits the knowledge of source posi- 
tions to extract as much information as possible on their 
fluxes. 

Since point sources, as any other foreground emission, 
are a nuisance for CMB experiments, these are designed 
to keep them as low as possible. Thus, if we are interested 
in point sources, it is of great importance to develop and 
apply specific and highly efficient detection algorithms. 
It has been sh own that wavelet techniques perform ve ry 
well t his task (ICavon et al.1 120001: IVielva et a!ll2001allE 
2003 ; [Gonzale z- Nuevo et al.l l2006t lLopez-Caniego et al.l 



2006). 



In a recent paper (jGonzalez-Nuevo et al.|[2006h some of 
us have discussed a natural generalization of the (circu- 
lar) Mexican Hat wavelet on R 2 , obtained by iteratively 
applying the Laplacian operator to the Gaussian func- 
tion, which we called the Mexican Hat Wavelet Fam- 
ily (MHWF). We demonstrated that the MHWF per- 
forms b etter than the stan dard Mexican Hat Wavelet 
(MHW, ICavon et al.l [2000f) for the detection of EPS 
in CMB anisotropy maps . In a subsequent work 
(|L6pez-Caniego et all I2006Q the MHWF has been ap- 
plied to the detection of compact extragalactic sources 
in simulated CMB maps. In that work it was shown, in 
particular, that the second member of the MHWF, called 
the MHW2, at its optimal scale compares very well with 
the standard Matched Filter (MF); its performances are 
very similar to those of the MF and it is much easier to 




Fig. 1. — S ky coverage of 5 GH z surveys in eq uatorial coor- 
dinat es: GB6 (IGregorv et al.H1996ri (bluel. PMNE (IGriffith et al.l 
IT995T) Cdark greenl PMNS flWright et all 11991) (red), PMNT 
(IGriffith et~aT][T994h (light blue) and PMNZ HTWrignt et al.lll996h 
(magenta). The white regions are 'holes' in these surveys, that 
have been covered exploiting the NVSS and the SUMSS. 



implement and use 1 . We have therefore decided to apply 
the MHW2 technique to estimate the flux of the EPS in 
the 3-year WMAP maps for a complete sample of sources 
selected at lower frequencies. 

The outline of the paper is as follows. In Section [21 
we introduce the low-frequency selected sample we have 
used. In Section [3] we present and discuss the tools for 
detecting the EPS and for estimating their fluxes. In 
Section 3] we present and describe our final catalogue. In 
Section[5]we compare our catalogue with the WMAP one. 
Finally, in Section[51 we summarize our main conclusions. 

2. THE 5 GHZ CATALOGUE 

A summary of the multi-steradian surveys that we 
have used is given in Table [TJ The highest frequency for 
which an almost complete sky coverage has been achieved 
is ~ 5 GHz, thanks to the combined 4.85 GHz GB6 and 
PMN surveys with an angular resolution of 3.5' and 4.2', 
respectively, and a flux limit ranging from 18 to 72 mJy 
The sky coverage of these surveys is illustrated in Fig. [TJ 
Deeper and higher resolution surveys have been carried 
out at 1.4 (NVSS, Condon et al. 1998; FIRST, Becker et 
al. 1995) and 0.843 GHz (SUMSS, Mauch et al. 2003); 
altogether these surveys cover the full sky. 

As extensivel y discussed by many authors in the 
recen t past dToffolatti et all 119991: iDe Zotti et all 
1991 iMason et al l 120031 : iBennett et all [2001 : 
Henkel & Partridge! [20051 "flat-spectrum" AGNs 
and QSOs, i.e. sources showing a spectral index a ~ 
(S(u) oc v~~ a ), are expected to be the dominant source 
population in the range 30-100 GHz, whereas other 
classes of sources, and in particular the steep-spectrum 
sources increasingly dominating with decreasing fre- 



1 LoDcz-Canicgo ct al. (2006) compared the two techniques, MF 
vs. MHW2, by exploiting realistic simulations of CMB anisotropics 
and of the Galactic and extragalactic foregrounds at the nine fre- 
quencies, between 30 and 857 GHz, of the ESA Planck mission, 
considering the goal performances of the Planck Low and High 
Frequ ency Instruments, LFI and HFI (see Planck Collaboration, 
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quency, are only giving minor contributions to the 
number counts at W MAP frequencies and sensitivities 
(|De Zotti et alj 120061) . We therefore chose to adopt 
5 GHz as our reference frequency and used lower 
frequency surveys to fill the "holes" at 5 GHz. 

Altogether, the catalogues listed in Table Q] contain 
over 2 million sources, but we already know, from the 
analysis of the WMAP team, that for only a tiny fraction 
(~ 2 x 1CT 4 ) of them the WMAP data can provide useful 
information. Applying the MHW2 at the positions of 
all these sources would not only be extremely inefficient, 
but plainly unpractical because of the huge CPU time 
and disk storage requirements. Therefore, we decided 
to work with a complete sub-sample containing sources 
with S5GHZ > 500 mjy. This limiting flux corresponds to 
about 2-3 times the mean noise in the filtered images we 
will be dealing with (see §[3]). To fill the 5 GHz "holes" 
we have picked up NVSS or SUMSS (in the region not 
covered by the NVSS) sources brighter than 500 mjy at 
the survey frequency. In this way we obtained an all-sky 
sub-sample of 4050 objects, whose spatial distribution, 
in Galactic coordinates, is shown in Fig.[2j After having 
removed sources in the strip |6| < 5°, and in the LMC 
region (i.e. inside the circle of 5.5° radius centered at 
a = 5 h 23 m 34 s .7, S = -69°45'22", J2000; I = 280.47°, 
b = —32.89°) and the Galactic sources outside of these 
zones (Taurus A, Orion A & B, and the planetary nebula 
IC418/PMNJ0527-1241) we are left with 2491 objects 
making up our "Input Catalogue" (IC). 

A c ross-correlation of th e IC with the WMAP cata- 
logue (|Hinshaw et al.ll2006l ) with a search radius of 20.8', 
equal to the dispersion of the Gaussian approximation 
of the beam of the lowest resolution WMAP channel (23 
GHz), showed that 298 of the 323 WMAP sources have 
a counterpart in the IC. The other 25 WMAP sources 
(called missed sources) must be unusually weak at low 
frequency, either because have an inverted spectrum or 
are strongly variable and were caught in a bright phase 
by WMAP. They are thus interesting targets for further 
study, and we have investigated them too. As they have 
a different selection, these sources are listed separately 
from the others (Table 4) . 

3. METHODOLOGY 

The strategy used by the WM AP team to obtain th e 
catalogue of 323 point sources (|Hinshaw et al.l l2006f ). 
summarized in §[U used some approximations: 

- The candidate detections were selected as 5<r peaks, 
where a is the rms noise defined globally. However, 
the removal of the Galactic emissions by compo- 
nent separation methods is not perfect, and leaves 
non-uniform contributions to the noise. Clearly, 
the detection, identification and flux estimation of 
EPS can benefit from a local treatment of the back- 
ground. Furthermore, considering a global rms 
noise can lead to an underestimate of the error in 
regions of strong Galactic emission. 

- The intensity peaks were fitted to a Gaussian pro- 
file. However, if the real beam response function is 
non-Gaussian, a Gaussian fit can lead to systematic 
errors in the flux estimate. 

- Confusion due to other point sources that are close 



to the target one, albeit a less relevant effect, is an- 
other source of error that makes necessary to study 
the data locally. 

As mentioned above, we used, locally, the MHW2 that 
proved to be as efficient as the matched filter but easier 
to implement and more stable against local power spec- 
trum fluctuations (Gonzalez-Nuevo et al. 2006; Lopez- 
Caniego et al. 2006). The real symmetrized radial beam 
profiles given by the WMAP team have been used rather 
than their Gaussian approximations. The error on the 
source flux density is calculated locally as the rms fluctu- 
ations around the source. Finally, we also corrected the 
flux of all the sources for the Eddington bias, adopting a 
Bayesian approach (see § 13. 4| . 

3.1. From the sphere to flat patches: rotation & 
projection 

In order to avoid CPU and memory expensive iterative 
filtering in harmonic space we chose to work with small 
flat sky patches. For every source position in the IC 
and for each WMAP frequency we obtained a flat patch 
of approximately 14.6° x 14.6° size, by projecting the 
WMAP full-sky maps. The adopted pixel area is 6.87 x 
6.87 arcmin 2 (NSIDE=512), so that the patches are made 
of 128 x 128 pixels. The patch making goes as follows: 

- Given the source coordinates we ob tain the cor- 
respo nding pixel in the HEALPix (|G6rski et al.l 
120051 ) scheme. 

- The image is rotated in the ai m space so that 
the position of the point source is moved to the 
equatorial plane. This is done in order to mini- 
mize the distortion induced by the projection of the 
HEALPix non-square pixels into flat square pixels. 

- The pixels in the sphere in the vicinity of the centre 
are projected using the flat patch approximation 
and reconstructed into a 2D image in the plane. 

- The units of the images are converted from mK 
(WMAP units) to Jy/sr and finally to Jy, using 
the real symmetrized radial beam profiles to do the 
integrals over the beam area. 

The rotation/projection scheme described above intro- 
duces some distortions in the projected image. There are 
two different effects to be considered: a) the distortion 
introduced by the projection (from HEALPix pixels to 
flat pixels in the tangent plane); b) the distortion caused 
by the rotation (from HEALPix pixels to a^ m and then 
to HEALPix pixels again). We have studied these dis- 
tortions by simulating 840 sources having flux densities 
ranging from 500 mjy to 20 Jy with the real WMAP 
beam profile. First they are simulated without noise and 
then they are added on the combined WMAP maps, plac- 
ing them at different Galactic longitudes and latitudes. 
Applying to each of the simulated sources the same ro- 
tation/projection procedure as for real sources we have 
found that, as expected, the projection effects are small 
at the image center (near the tangent point) and grow 
towards the borders of the patch. Also distortion effects 
are small near the equator (where the HEALPix pixels 
are very close to squares) and grow towards the poles. 
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Fig. 2. — Sky distribution, in Galactic coordinates, of the 4050 objects with flux densities above our selection threshold. 



This is the main reason for performing the rotation in 
order to always have the point source at the equator. 
We find that the rotation in the ai m space from the ini- 
tial position of the source to the centre of the map (1=0, 
6=0) has a very small effect on the flux estimate (always 
< 2%), while the effect of the projection is totally negli- 
gible. 

3.2. Filters 

In pre vious works (iSanz et al.l 11999b ICavon et al.l 
20001: ISanz. Herranz. & Martinez-Gonzalez! l200lh 
Vielva et all l2001alfbl. l2003t iHerranz et alJ l2002allblR 
Gonzalez- Nuevo et al.l 120061 : lLopez-Caniego et al.l 2006) 
we have shown that application of an appropriate filter 
to an image helps a lot in removing large scale variations 
as well as most of the noise. As a result, the signal- 
to-noise ratio is increased, thus amplifying the point 
source signal. A brief summary of the fundamentals of 
linear filtering of two-dimensional (2D) images is given 
in Appendix In this work we have only taken into 
account two possible filters: the Matched Filter (MF) 
and the second member of the Mexican Hat Wavelet 
Family (MHW2). By applying these two filters to 
WMAP temperature maps we find that both of them 
give an average amplification of the EPS flux by a factor 
of ~ 3. This is a remarkable result. It means that the 
flux of a point source at the ~ 2a level in the original 
map can be enhanced to a > 5<r level in the filtered 
map. In the following sub-sections, we briefly sketch the 
main properties of the MF and of the MHW2. 

3.2.1. The Matched Filter (MF) 

The MF is a circularly-symmetric filter, ^>(x;R, b), 
such that the filtered map, w(R, b) (see Appendix A for 
the definition of the notation) , satisfies the following two 
conditions: (1) (w(Rq,0)) = s(0) = A, i.e. w(R, 0) is an 
unbiased estimator of the flux density of the source; (2) 
the variance of w(R, b) has a minimum on the scale Rq, 
i.e. it is an efficient estimator. In Fourier space the MF 
writes: 

where P(q) is the power spectrum of the background and 
r(q) is Fourier transform of the source profile (equal to 
the beam profile for point sources). The matched filter 
gives directly the maximum amplification of the source 



and it yields the best linear estimation of the flux, when 
used pro perly and under contr o lled c onditions. As men- 
tioned in lLopez-Caniego et alJ (|2006h . the practical im- 
plementation of the MF requires the estimation of the 
power spectrum P(q) directly from the data and this 
leads to a certain degradation of its performance. 

The WMAP team have done a global implementation 
of the MF on the sphere by taking into account the non- 
Gaussian profile of the beam (although the source fluxes 
are estimated fitting the source profiles with a Gaussian) . 
The resulting matched filter is given by eq. ([T]), where 
the flat limit quantities r(q), P(q) are replaced by their 
harmonic equivalents be, Ce- The use of the Ce of the 
whole sky to construct a MF filter that operates in the 
sphere is a good first approach to obtain a list of source 
candidates and to estimate their fluxes, but we do believe 
that it can be improved by operating locally. In this 
paper we use different filters for regions with different 
levels of Galactic contamination. 

3.2.2. The Mexican Hat Wavelet Family (MHWF): MHW2 

One example of wavelet that is particularly well suited 
for point source detection is the MHW2, a member of 
the M HWF first introduced by iGonzalez-Nuevo et al.l 
(|2006h . The MHW2 is obtained by applying twice 
the Laplacian operator to the Gaussian function. This 
wavelet filter operates locally and is capable of remov- 
ing simultaneously the large scale variations introduced 
by the Galactic foregrounds as well as the small scale 
structure of the noise. Note that the expression of the 
filter can be obtained analytically (while the MF depends 
on the P(q) that must be estimated numerically and lo- 
cally 2 ). Any member of this family can be described in 
Fourier space as 

i> n (q) cx q 2n e-<?l\ (3) 

2 Some controversy has arisen lately on whether the dependence 
of the MF upon P(q) is a problem from the practical point of view, 
or not. The main argument supporting the opinion that there is no 
problem at all goes as follows: the angular power spectrum of the 
background signal is determined by the CMB, whose power spec- 
trum is fairly well known, and by the instrumental noise, whose 
statistics is perfectly known. This is not quite true since the back- 
ground signal includes Galactic emission (or its residual after com- 
ponent separation), which shows strong variations from one point 
of the sky to another, and unresolved point sources. Thus the 
angular power spectrum is not perfectly known, at least locally. 
It must be estimated from the data with some error that will in- 
evitably propagate to the fil ter. For a more detailed discussion, see 
Lopez-Caniego et al. (2006). 
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Fig. 3. — Example of a bright source filtered with the MHW2 (left 
in the case of the MF. 

The expression in real space for these wavelets is 

1> n (x) oc A>(i), (4) 

where tp is the 2D Gaussian <p(x) = (l/2ir)e- x ' 2 . We 
remark here that the first member of the family, tpi(x), 
is the usual Mexican Hat Wavelet (MHW), that has been 
exploited for point source detection with excellent results 
(jCavon et al.ll2000l ; IVielva et al.ll2001all2003h . Note that 
we call MHWn the member of the MHWF with index n. 

As already mentioned, in this work we filter our pro- 
jected WMAP sky patches with the second member 
of the family, the MHW2. As in lLopez-Caniego et al.l 
we will do a qualitative comparison with the re- 
sults obtained with the MF, implemented to be used lo- 
cally in flat patches (at variance with the global MF used 
by the WMAP team). 

3.3. Position, flux and error estimation 

We want to obtain an estimate of the flux density, with 
its error, of the IC sources at the center of each filtered 
image. Point sources appear in the image with a profile 
identical to the beam profile. For example, if the beams 
were Gaussian, the ratio between the beam and the pixel 
area, 2ir(R s ) 2 /Ll, where R s = FWHM/(2V21og2) is the 
beam width and L p is the pixel side, would allow us to 
convert the flux in the pixel where the source is located 
into the source flux. 

In our case, the beams are not Gaussian and we will 
calculate this relationship integrating over the real beam 
profile for each channel. In carrying out the calculation 
we have to take into account that we work with HEALPix 
coordinates, at the WMAP resolution. Although the im- 
age is centered on the source position, after the projec- 
tion to the flat patch the source does not always lay in 
the central pixel, but may end up in an adjacent one. 
Thus, to estimate its flux we make reference not to the 
intensity in the central pixel but to that of the brightest 
pixel close to the center of the filtered image (however in 
most cases the brightest pixel coincides with the central 
one). 

At first glance, this may seem a very crude estimator, 
but it turns out that flux estimation through linear fil- 
tering is almost optimal in many circumstances. Let us 
explain how this estimator works. After filtering, the in- 
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panel) and the MF (right panel). Strong ringing effects can be seen 

tensity of the brightest pixel can be written as a weighted 
sum of the intensities in the surrounding pixels, 

w(x ) = ^ip(x k - x )Ip(x k ) (5) 

k 

where xq is the position of the considered source, I p {xk) 
is the intensity of the pixel of the unfiltered image located 
at the position Xk and ip is the kernel of the filter. It can 
be shown that if tp = tJjmf (that is, for the matched 
filter) w(xo) is the best possible linear estimator (statis- 
tically unbiased and of maximum efficiency) of the flux of 
the source. In particular, for the case of Gaussian noise 
w(xo) is the maximum l ikelihood estimator of the flu x of 
the source. As shown in lLopez-Caniego et al.l (|2006l) the 
flux estimation when tp = ip2 (that is, with the MHW2) 
is comparable to that obtained with the MF. 

The method adopted here differs from the one used 
by the WMAP team. They have used the MF to detect 
point sources above the 5tr level in the filtered images 
(full-sky maps, in their case), but not to estimate their 
fluxes. These are derived by fitting, in the unfiltered 
image, the pixel intensities around the point source to 
a Gaussian profile plus a plane baseline. However, as 
already noted, the profile of the source matches that of 
the beam, and is therefore non-Gaussian. In this paper 
the true source profiles are used at every step, up to the 
final flux unit conversions. 

For each source we calculate the dispersion a as the 
square root of the variance of the background pixels that 
are close to the target source but are not "contaminated" 
by it. This number can be easily inflated in the pres- 
ence of contamination by nearby sources or large scale 
structures that, in some cases, cannot be removed com- 
pletely by filtering. In order to avoid this, we first select 
a shell of pixels around the source, with an inner radius 
equal to the FWHM of the beam. This guarantees that 
the source flux has decreased well below the background 
level. Subsequently, we choose an outer radius encom- 
passing a sufficiently large number of pixels (~ 3000) to 
give an accurate estimate of a. As a final step, we di- 
vide this shell in sectors - 12, normally - and calculate 
their mean and dispersion. Strong contamination, due 
to other sources present in the annulus, shows up as a 
significant difference between the mean and the median; 
whenever this difference exceeds two times the dispersion 
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Fig. 4.-5 GHz vs 23 GHz WMAP flux densities of the 25 
sources missed by our low-frequency selection. The source with 
zero 23 GHz flux was not detected at 23 GHz, but at least at one 
higher frequency. Encircled dots are sources that we recover at 
> 5a. 

we excluded the pixels in such sector from the calculation 
of the variance. 

The mean values of a, the rms error on the EPS fluxes, 
turn out to be: 198, 231, 222, 255 and 399 mjy at 23, 
33, 41, 61 and 94 GHz, respectively. The variation of 
a with the WMAP channel is determined by the beam 
shapes, the spectral behavior of the foreground emission, 
the instrumental noise, etc.. It may be noted that these 
uncertainties are typically 2 or 3 times higher than the 
uncertainties quoted by Hinshaw et al. (2006), which are 
probably underestimated, as confirmed by the fact that 
their source counts show clear signs of incompleteness at 
flux densities well above 5 times their typical errors. 

3.4. Bayesian correction to the fluxes 

Flux-limited surveys are affec ted by the well-known 
Eddington bias (Eddington 1940) that leads to an overes- 
timate of the flux of faint sources. Hogg & Turner (1998) 
have shown how to correct for this effect if the underly- 
ing distribution of source fluxes is known. Unfortunately, 
our knowledge of this distribution in the frequency range 
covered by WMAP is rather poor. However, as shown in 
the Appendix |BJ given a set of observed fluxes and the 
associated values of the rms noise it is possible to simul- 
taneously estimate the slope of the flux distribution (i.e., 
of the differential number counts) and to obtain an unbi- 
ased estimate of the source fluxes. We have applied this 
method, based on th e Bayesian approach introduced by 
iHerranz et al.l (|2006l ). to correct for the Eddington bias. 
The estimated slopes of the differential number counts 
are 2.11, 2.34, 2.16, 2.14, and 2.16 for the 23, 33, 41, 61, 
and 94 GHz channels, respectively, in very good agree- 
ment with the results of the ATCA 18 GHz survey (Ricci 
et al. 2004: 2.2 ±0.2), of the 9C survey at 15 GHz (Wal- 
dram et al. 2003: 2.15), and of the 33 GHz VSA survey 
(Cleary et al. 2005: 2.34±g:||). 

4. THE NEW EXTRAGALACTIC WMAP POINT SOURCE 
(NEWPS) CATALOGUE 

4.1. Comparison between MHW2 and MF 



We have carried out the detection/flux estimation pro- 
cess using the two different filters, MHW2 and MF, pre- 
viously discussed. We have found that filtering with the 
MF introduces ringing effects around the target EPS 
in at least 15% of the images analyzed (see Figure [3|). 
These effects are stronger for the brightest sources. Rings 
appear due to strong oscillations in the shape of the 
matched filter which are determined by sharp features in 
the angular power spectrum of temperature anisotropies 
in the sky patch. These sharp features are likely to ap- 
pear in regions showing a high background signal and/or 
in regions where point sources constitute a major com- 
ponent of the total intensity of the image. If the num- 
ber of images were small, we could check the images 
one-by-one visually and study separately the anomalous 
cases, but this procedure is impractical in the present 
context. With a blind approach we get a number of spuri- 
ous sources, mainly due to positive interferences between 
rings. In the case of our non-blind approach rings may 
also contaminate surrounding pixels, thus affecting the 
error estimates and the flux estimates of the surround- 
ing sources. On the contrary, we found that the MHW2 
filter does not introduce ringing effects in any of the con- 
sidered images. 

The number of source detections above a certain no~ 
threshold obviously depends on the correct estimation of 
the background noise level. Ringing effects affect nega- 
tively such estimate. We find that, on average, by fil- 
tering with the MHW2 we correctly identify ~ 7% more 
sources than with the MF 3 . For this reason we decided 
to build the final catalogue using the detect ions /flux es- 
timates obtained with the MHW2. 

4.2. NEWPS 3a 

To define our first SubCatalogue (SC), NEWPS 3<T , of 
sources at > 3<r we need to take into account that the 
IC was generated from surveys with angular resolution 
< 4'. 2, i.e. much higher than the WMAP resolution, so 
that we may have multiple IC sources within one WMAP 
beam. In some cases the WMAP signal is entirely due 
to the brightest source within the beam, whose flux, at- 
tenuated by the response function, accounts also for the 
signal detected in the direction of nearby sources. The 
latter are therefore false detections that we have removed 
from the sample. A search within circles of radius equal 
to 2 WMAP beam sizes (beam sizc=FWHM/2 y/2 log(2)) 
centered on the position of each source, starting from the 
brightest ones, has yielded numbers of false detections 
decreasing from 354 at 23 GHz to 104 at 94 GHz. 

The contamination of very bright sources can actually 
extend beyond 2 WMAP beam sizes. We have therefore 
repeated the above search up to 4 beam sizes, and re- 
moved those sources for which the contamination by the 
bright source accounts for more than 50% of the detected 
signal. The number of such cases oscillate around 7 at 
all channels. In 13% of the cases, the data are consistent 
with more than one IC source contributing significantly 
to the WMAP flux. 

After having removed the 4 |6| > 5° sources known 
to be Galactic (Taurus A, Orion A & B, and IC418), 
we are left with 760(7384-22 missed sources, WMAP 

3 This percentage varies depending on the observation frequency, 
from ~ 14% to ~ 4%. 
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Fig. 5. — Comparison of N EW PS^a with WMAP flux densities at 23 and 61 GHz: Blue circles show our uncorrected fluxes and red 
asterisks the corrected ones. 




Fig. 6. — Number counts multiplied by S 5 / 2 . Red asterisks: 
counts from the N EW PS§ a catalogue at 33 GH z (see text for 
more details). Black diamonds: WMAP counts (Hinshawetal 
2006]) . Black dots: ATCA 18 GHz pilot survey counts IIRicci et al 
2004 ). The paralle logram is from the DASI experiment at 31 GHz 
pCovac et all 12001 ). The solid c urve shows, for comp arison, the 
counts predicted by the model by Do Zotti et al. (2005). 



sources that are not included in our IC), 565(547+18), 
536(518+18), 366(360+5) and 103(101+2) EPS with 
signal-to-noise ratios > 3 in the 23, 33, 41, 61, and 94 
GHz channels, respectively (see Table[5]). Finally, a cross 
correlation of the catalogues at the different frequencies 
shows that we have 933 (908+25) different EPS with 
signal-to-noise ratios > 3 in at least one WMAP channel, 
which constitutes our N EW P S^ a Catalogue. 4 

4.3. NEWEST : A robust subsample 

Of the 381 sources that we detected at > 5a 
(NEWPS 5cT Catalogue in Table 4, see Appendix C, 
plus 12 EPS detected at > 5a listed in Table 5) only 
283 are listed in the WMAP 3-yr catalogue (see Ta- 
ble [3] for a better understanding of the detection statis- 
tics). As expected, the a priori knowledge of source po- 
sitions has allowed us to significantly increase the detec- 
tion efficiency. Additionally, in our approach, 39 (26+13) 
WMAP sources present a signal-to-noise ratio < 5 (as 
shown in Table 3) 5 . 

4 All the catalogues can be obtained from 
max.ifca.unican.es / caniego /NEWPS 

5 However, 3 of the 26 WMAP sources with low frequency coun- 
terparts included in the IC are not 5<r detections in the original 
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Regarding the reliability of the detected sources near 
the 5a threshold, we would expect that most of them, if 
not all, are actually real 5a sources. Since we have fol- 
lowed a non-blind approach we know that there is a real 
source in every position we have considered, but it is pos- 
sible that a very weak source appears above this thresh- 
old because it has been contaminated, e.g. by Galactic 
foregrounds, CMB, etc. In order to check this possibility 
we have created a catalogue of false positions and we have 
repeated the process of rotation, projection, filtering of 
each patch and flux estimation exactly in the same way 
as we have done for detecting the sources in NEWPS$ a . 
Then we have counted the number of > 5a detections. 
This number gives an idea of how likely is that the fore- 
grounds produce a spurious > 5a detection and therefore 
how reliable is our catalog for that threshold. 

This catalogue of false positions is obtained by shift- 
ing 10 degrees in the longitude coordinate the position 
in the sky of each source outside |6| > 5° in the IC. 
We have performed this test only at 23 GHz, because at 
this frequency we have the highest number of detections 
and the lowest mean value of rms noise. As a result, we 
have found 3% of spurious > 5a detections in the anal- 
ysed patches, which is in agreeme nt wit h the numbers 
discus sed by iBennett et al.l (|2003l ) and iHinshaw et al.1 

(pool . 

In addition, as already mentioned, there are 25 WMAP 
sources left out by our low-frequency selection, but in- 
cluded in our analysis (see Table 4). Only 12 of them 
are detected at > 5tr by our approach, and only 3 have 
23 GHz flux densities above the estimated completeness 
limit of 1.1 Jy (see Fig. [6]). Of our 12 sources which 
are > 5a detections, 10 have low-frequency flux densities 
> 270 mJy (8 of them are above 340 mJy) and may well 
be variable sources, that happened to be in a particu- 
larly 'high' phase at the time of WMAP observations. 
The source with Ssghz = 120 mJy has an inverted spec- 
trum (i.e. a spectrum rising with frequency) and the last 
one, with Ssghz = 38 mJy, is a strong candidate to be a 
spurious source. 

Finally, by exploiting the distribution of spectral in- 
dexes a| 3 , we could also tentatively estimate how many 
more 5a detections we missed at 23 GHz because of the 
adopted 5 GHz flux limit S$ghz > 500 mJy. From 
the observationally determined number counts at 5 GHz 
(and the associated uncertainties) we estimate that there 
should be ^2100-2500 sources in the flux range 300 < 
S^GHz < 500 mJy, in the same sky area in which we find 
2491 sources with S > 500 mJy. For being detectable at 
23 GHz, these sources should show an inverted spectrum 
with slope a < —0.63. The distribution of the spectral 
indexes a^ 3 of sources in the IC shows that only < 1% of 
them show spectra so inverted. This implies that we may 
expect < 20-25 sources with 300 < S^ghz < 500 to be 
detectable at > 5a at 23 GHz. The number of expected 
detections decreases to < 5 for 100 < S^ghz < 300, and 
is essentially zero at still fainter 5 GHz fluxes. These esti- 
mates are comparable with the number (11) of sources at 
23 GHz (see Table 3) detected by WMAP but missed by 
our selection criterion. Note that, at higher frequencies, 
the numbers of missed detections, which is more difficult 

WMAP catalogue. These sources are: GB6 J1228+1124, GB6 
J1231+1344 and GB6 J1439+4958. 



to estimate due to the higher errors in flux, decreases 
in parallel with the decrease of the number of detections 
(see Table 2 and Section 5.1). 

5. RESULTS AND DISCUSSION 
5.1. NEWPS 5a versus WMAP sources 

The WMAP 3-year catalogue (Hinshaw et al. 2006) 
lists 314, 292, 280, 154, and 29 > 5a detections at 23, 
33, 41, 61, and 94 GHz, respectively, to be compared with 
350 (339+11 missed sources), 224 (219+5), 218 (215+3), 
136 (135+1) and 22 sources in the NEWPS 5a catalogue 
(see Table [2]). Although the blind WMAP technique 
yielded a lower number of detections at 23 GHz and a 
lower total number, it was apparently more successful in 
retrieving sources at > 5a in different channels. This 
may be, to some extent, related to their lower error es- 
timates. However, note that the same sources can be 
detected in different channels. Therefore, the total num- 
ber is the intersection of the channels, giving a global a 
reduced number for the WMAP case (323 vs. 381). 

In Fig. [5] we compare our flux estimates for sources 
in the NEWPS 3r7 catalogue with the WMAP ones at 2 
frequencies (23 and 61 GHz). We have plotted the cor- 
rected fluxes (red asterisks) as well as the uncorrected 
ones (blue circles) to show the effect of the Bayesian cor- 
rection. At 23 GHz the agreement is generally good and 
the correction makes no difference for fluxes > 1.5 Jy, 
due to the fact that almost all the sources are detected 
at high signal to noise ratios. 

The most striking difference is found for Fornax A 
(PMNJ0321-3658), represented by the isolated point on 
the top right-hand corner of the 23 GHz channel and at 
the center of the upper part of the 61 GHz panel. Our 
procedure yields 23 GHz and 61 GHz flux densities of 
9.1 + 0.2 and 1.6 + 0.2 Jy, to be compared with 18.5 + 3.6 
and 8.3 ± 2.1 Jy, respectively, given in the WMAP cat- 
alogue (fluxes obtained by aperture photometry). The 
problem here is that Fornax A has a relatively weak core 
and 2 big lobes extending, altogether, over ~ 1°, so that 
it cannot be treated as a point source, even at the rela- 
tively low resolution of WMAP. As both flux estimates 
are not reliable we have applied the NRAO Astronomical 
Image Processing System (AIPS) software to the unfil- 
tered patches, fitting the source to a gaussian. The values 
obtained by this method are 11.49 ± 0.63 and 0.92 ± 0.72 
Jy, with a best fit major(minor) axis of 1.18°(0.81°). 
Resolution effects may be worrisome also for Centau- 
rus A (PMNJ1325-4257), the largest radio source in the 
sky at low frequencies, which is however not present in 
the WMAP catalog. Therefore for this source too we 
have checked our flux estimates with those obtained from 
AIPS, which gave 48.12 + 0.04 and 20.31 + 0.72 Jy, for 23 
GHz and 61 GHz respectively, and a best fit major (mi- 
nor) axis of 1.13°(0.84°) for 23 GHz, in good agreement 
with our results. 

The second discrepant point (circle not far from the 
panel center at 23 GHz) is PMN J0428-3756 for which 
we get a 23 GHz flux of 1.12 Jy, while the WMAP flux is 
2.20 Jy. This source is an ATCA flux calibrator. Its light 
curve 6 at 20 GHz shows its flux increasing from ~ 1 Jy 
in 2002 to a maximum of 2.04 Jy in 2004, and steadily 

6 www.narrabri.atnf.csiro.au/cgi-bin/Calibrators/calfhis.cgi? 
source=0426-380&band= 1 2mm 
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TABLE 1 

Summary of the large-area surveys of point sources used in this work. 



Frequency Catalogue Su m (m3y) 



DEC range 



Angular resolution 



4.85 GHz 



GB6 
PMNE 
PMNT 
PMNZ 
PMNS 



18 
40 
42 
72 
20 




-9.5 
-29 
-37 
-87.5 



+75 
+10 
-9.5 

-29 

-37 



3.5' 
4.2' 
4.2' 
4.2' 
4.2' 



Gregory ct al 
Griffit h et al. 
Griffit h et al. 
Wrigh t et al. 
Wright et al. 




Condon ct al (1998)~ 
Mauch ct al. (2003) 



1.4 GHz 



NVSS 



2.5 



-40 



+90 



45" 
45" coscc 1 5 1 
45" coscc 1 8 1 



0.843 GHz 



SUMSS 



18 
8 



-50 
-90 



-30 
-50 



TABLE 2 

Summary of the results. 





23 GHz 


33 GHz 


41 GHz 


61 GHz 


94 GHz 


Total 


< rms > 


198 


231 


222 


255 


399 




a 


2.11 


2.34 


2.16 


2.14 


2.16 




NEWPS Za 


760 


565 


536 


366 


103 


933 




(738+22) 


(547+18) 


(518+18) 


(360+6) 


(101+2) 


(907+25) 


NEWPS 5 „ 


350 


224 


218 


136 


22 


381 




(339+11) 


(219+5) 


(215+3) 


(135+1) 


(22+0) 


(369+12) 


WMAP 5a 


314 


292 


280 


154 


29 


323 



Note. — The "rms" of the patches is in units of mjy. The slope a corresponds to 
the estimated slope of the flux distribution (i.e. of the differential number counts, dN/dS 
oc S~ a ). For the N EW PS$<j and N EW PS^^ we show in parenthesis the number of 
detections coming from the initial catalogue "IC" and the number of detections among the 
25 WMAP objects not present in IC (missed sources, see Table 4). 



TABLE 3 
Detection statistics 





(All v) 


(All v) 


(23 GHz) 


(23 GHz) 


SubCatalogue (SC): 


3a 


5a 


5a 


> Si im i-Uy 


Number of WMAP EPS 6 SC 


297 


271 


258 


186 


Number of WMAP EPS g SC 




26 


39 


111 


Number of Missed EPS 6 SC 


25 


12 


11 


3 


Number of Missed EPS g SC 




13 


11 


22 


Dropped WMAP EPS 


1 


1 


1 


1 


New EPS in WMAP) 


611 


98 


81 


13 



Note. — Distribution of the detected sources in different subsamples. The de- 
scription has been made for 4 classes of EPS: WMAP EPS (WMAP sources that 
appeared in our IC), Missed EPS (WMAP EPS NOT included in our IC), Dropped 
EPS (the planetary nebula IC418 (PMNJ0527-1241)) and New EPS (detected EPS 
NOT included in the WMAP EPS catalogue). The total number of sources in each 
SubCatalogue, detected by the MHW2 filter, are given by the sum of the entries in 
lines (l) + (3) + (6), in agreement with the numbers in Table 2. 

decreasing afterwards back to ~ Uy in 2006. 5.2. Source number counts 

In Figure [6] we compare the number counts derived 
from the N EW PSz a catalogue at 33 GHz (red asterisks 
with Poisson error bars) with other sets of observational 
data, specified i n the caption, a nd w ith the prediction 
of the model by iDe Zotti et alj (|2005|) . The agreement 
is clearly good. The completeness flux limit of our cat- 
alogue is around ~1.1 Jy, as expected from the average 
value of the flux error at this frequency, ~ 230 mjy. 

Hinshaw et al. (2006) estimated the contribution 
of point sources below the detection threshold to the 
anisotropy power spectrum at 40.7 GHz to be A — 
0.017 + 0.002 /jK 2 . Again, this is in good agreement with 
the De Zotti et al. (2005) model, yielding, at this fre- 
quency, A = 0.018 nK 2 . 



The asterisks in the lower left-hand part of each panel 
correspond to objects that we detect at < 5cr at 23 GHz, 
whose flux density is substantially decreased by the cor- 
rection for the Eddington bias. 

At 61 GHz WMAP flux densities are systematically 
brighter than ours by ~ 10%. This may be due to the use 
by the WMA P team of Gaus s ian so urce profiles while, as 
discussed bv iHinshaw et al.1 (|2006h . beams increasingly 
deviate from Gaussianity with increasing frequency. On 
the other hand, we have taken into account real sym- 
metrized beam profiles. At this frequency many sources 
fainter than ~ 2 Jy are detected by our procedure at < 5<j 
and have substantial corrections for the Eddington bias, 
moving them to the left of the diagram. 
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6. CONCLUSIONS 

We have used the MHW2 filter (Gonzalez-Nuevo et al. 
2006) to obtain estimates of (or upper limits on) the flux 
densities at the WMAP frequencies of a complete all-sky 
sample of 2491 sources at |6| > 5°, brighter than 500 mJy 
at 5 GHz (or at 1.4 or 0.84 GHz, in regions not covered 
by 5 GHz surveys but covered by either the NVSS or the 
SUMSS). We have shown that the MHW2 filter has an 
efficiency very similar to the MF and is much easier to 
use. 

We have obtained flux density estimates for the 933 
sources detected at > 3er, our NEWPS^a Catalogue 
(see Section 4.2). Of the 381, presumably extragalac- 
tic, sources detected at > 5cr, 369 of which constitute 
our NEWPS 5a Catalogue (Section 4.3), 98 (i.e. 26%) 
are "new" , in the sense that they are not present in the 
WMAP catalogue; 43 of them are above the estimated 
completeness limit of the WMAP survey: ~ 1.1 Jy at 23 
GHz. This illustrates how the prior knowledge of source 
positions can help their flux measurements. 

On the other hand, 39 (23+13) WMAP extragalactic 
sources with low-frequency counterparts included in our 
sample were detected by us at < 5o\ This is probably 
due to the fact that our error estimates exceed substan- 
tially those given by the WMAP team, particularly at 
23 GHz, where we have the highest detection rate. In 
fact, the WMAP errors do not correspond to the rms 
fluctuations in the source neighborhood but to the un- 
certainties in the amplitude of the Gaussian fit. As a 
consequence, the WMAP catalogue starts being incom- 
plete at flux densities more than 2 times higher than 5 
times their typical formal errors. 

Our flux density estimates for sources detected at > 5cr 
are generally in very good agreement with the WMAP 
ones at 23 GHz. At higher frequencies WMAP fluxes 
tend to be slightly but systematically higher than ours, 
probably due to having ignored the deviations, increasing 
with frequency, of the point spread function from a Gaus- 
sian shape. Our estimates use the real beam shape at ev- 
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ery frequency. For only one source we have a strong dis- 
crepancy with WMAP. Such source, Fornax A, is known 
to have powerful lobes extending over ~ 1°, and is there- 
fore resolved by WMAP. Thus, the point source assump- 
tion on which both WMAP and our flux estimates rely, 
is clearly not appropriate and yields unreliable results. 
A smaller, but still significant, discrepancy, is found for 
PMN J0428-3756. 

We have also worked out and applied a method to 
correct flux estimates for the Eddington bias, without 
the need of an a-priori knowledge of the slope of source 
counts below the detection limit. 

Our selection criterion leaves out 25 WMAP sources, 
only 12 of which however turn out to be > 5<r detec- 
tions after our analysis, and only 3 have 23 GHz fluxes 
> 1 . 1 Jy, the estimated completeness limit of the survey. 
Thus, our approach has proven to be competitive with, 
and complementary to the blind technique adopted by 
the WMAP team. In fact we missed 3 sources brighter 
than the estimated completeness limit (S*23GHz = 1-1 Jy), 
but we detected 42 new ones. 

On the whole, 26% of sources we have detected at > 5cr 
are not present in the WMAP catalogue. On the other 
hand, the efficiency of the process is low. Only 381 of the 
2491 sources in our input sample were detected at > 5a 
in at least one WMAP channel. 

In a forthcoming paper, we will exploit our catalogue 
to investigate the high-frequency properties of sources 
selected at low frequencies. 
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APPENDIX 

LINEAR FILTERS IN TWO DIMENSIONS 

The approach adopted in this paper assumes that in the center of each image (sky patch) there is a point source 
with unknown flux density, A. As usual we describe the source as 

s{x)=At(x), (Al) 

where t(x) is the source profile. We will assume circular symmetry, so that t(x) = t(x), x = \x\. For point sources 
the profile is equal to the beam response function of the detector. For a circular Gaussian beam we have 

T (x) =e- x2 / 2R ", (A2) 
where R s is the angular radius of the beam response function. In the Fourier space the beam profile of eq. (|A2|) is 

r(q) = e-?<l\ (A3) 

where q = \q\. Actually the WMAP beams are not Gaussian and we will use the real symmetrized radial beam profiles 
for the different WMAP channels to construct our filters. 
Let us consider a 2D-filter ^>(x; R, b), where R and b define a scaling and a translation respectively. Then 

9&R,$ = ^^±^y (A4) 

If we filter our field f(x) with ^(x; R, b), the filtered map will be 

w{R,b)= J dxf(x)f{x;R,b). (A5) 
The filter normalization preserves the amplitude of the source at its position (b = 0) after filtering: 

dxT(x)V(x;R,0) = 1. (A6) 
The n-th moment of the filtered map is defined as 



(7 



2 



2tt / dq q 1+2n P(q)lP 2 (q). (A7) 



where P(q) is the power spectrum of the unfiltered map and tpil) is the Fourier transform of ij)(x): 

/>oo 

ip(q) = / dxxJ Q (qx)ilj(x). (A8) 



Here q is the wave number, q = \q\ and Jo is the Besscl function of the first kind. The zeroth-order moment, o~o, which 
is the dispersion of the filtered map, will be denoted, for simplicity, as a. It is a function of the filter scaling R, of 
the source profile r and of the power spectrum P(q) of the unfiltered image. Then, for a given image and a given 
source profile it is possible to optimize the ratio A / a 'just by modifying the scaling R. The scale Rg at which A / a is 
maximum is called optimal scale. It has been shown (|Vielva et al.ll2001al 120031 : lLopez-Caniego et alJl2004L l2005| ) that, 
working at this scale at which a filter maximizes the flux of a compact source with respect to the average surrounding 
background, the filter performance in detecting point sources and estimating their flux is maximized. 

DERIVATION OF THE BAYESIAN CORRECTION FORMULAE 
The Eddington bias 

The distribution, normalized to unity, of true fluxes, S, of extragalactic sources, is usually well described by a power 
law 

P(S\q) = kS-( 1+q \ S>S m (Bl) 
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where the normalization k is 

k = qS m \ q>0. (B2) 

The observed fluxes are contaminated by noise. Let S° = {5°, . . . , S^} be the observed fluxed of the N galaxies 
detected above a given flux threshold and S = {Si, . . . , Sn} their true fluxes. In the case of Gaussian noise we have 

N , ™ r~f \ 2~1 



P(S°|S,«) oc exp 



4e 



S? Si 



(B3) 



where Ui is the rms noise for the ith source. Since we select only sources above a certain threshold, a selection effect 
appears: sources whose intrinsic flux is lower than the detection threshold may be detected due to positive noise 
fluctuations and, on the other hand, sources that should be detected because their flux is greater than the detection 
threshold may be missed due to negative fluctuations. Since according to eq. ()B1|) there are more faint than brigh t 
sources we end up with an excess of sources whose flux is overestimated. This is the Eddington bias (|Eddingtoni ri940) . 

Bayesian slope determination and flux correction 
^From the Bayes' theorem we have 

P(S°|S,g) P(S,q) 



p(s, q s°) = 



P(S°) 

P(S°|S,<?) P(S\q)P(q). 



(B4) 



We will assume that P(q) is uniform, that the source fluxes take values drawn at random from the distribution of 
eq. (IBip . and that the flux of any given source is independent of the fluxes of the other sources. Therefore 

P(S\q)=l[P(S i \q) = q N S^ [USA . (B5) 



i=i 



Therefore, using eqs. (|B3p and (|B4[) we have 

P(S,g|S°)cxexp 



2 ^ 

i=l 

i N s»« (f[s l 



S? Si 



-(9+1) 



(B6) 



If the slope q is known, the maximum likelihood estimator of the fluxes of the sources is easily calculated (Hogg 
Turner 1998): 



CO 

2 




4(1 + 9) 



> 2^1 + q, 



(B7) 



where = S°/ai is the signal to noise ratio of the source. Conversely, if the true intrinsic fluxes of the sources S were 
known and the slope q unknown, the maximum likelihood estimator for the value of q would be 

-l 



if, (Si 



i=l 



(B8) 



Unfortunately, in many cases, and particular in the case of WMAP surveys, neither q nor S are kn own a priori. The n 
it is necessary to solve simultaneously for the two unknowns. A way to do this is to introduce eq. (|B7j) into eq. (1B8|) . 
which gives the implicit equation 

i i f r, fs? \ r i 



1 + A, 



(B9) 



where 



4(1 + q) 



•l- 



(B10) 



Equation (|B9|) can be solved numerically if the minimum signal to noise ratio of the galaxies considered satisfies the 
condition > 4(1 + q). Once q is estimated, the fluxes S can be estimated using eq. (|B7[) . 
The asymptotic limits of the estimators, valid in the high signal to noise regime, that is, for r m 3> q, are: 

-l 

S° 

Co 
•J I, 



i=l 



Si — Sa 



1 



1 + 



(Bll) 
(B12) 



WMAP extragalactic point sources 
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TABLE 4 
Missed EPS 





RA 
h 




Dec 
deg 




K(a) 

Jy(Jy) 




Ka(cr) 

Jy(Jy) 




Q» 
Jy(Jy) 




V(cr) 

Jy(Jy) 




W(o-) 

Jy(Jy) 


Jy 


WMAP 


low freq. id. 


o 


434 


-35 


197 


1 


19 


(0.15) 





96 


(0.21) 


1 


25 


(0.20) 


1 


06 


(0.26) 





82 


(0.38) 


0.12 


6 


PMNJ0026-3512 


o 


493 


5 


922 


1 


09 


(0.16) 


1 


43 


(0.22) 





89 


(0.21) 





74 


(0.26) 


1 


22 


(0.45) 


0.35 


7 


PMNJ0029+0554 


5 


229 


-20 


268 





71 


(0.17) 


o 


51 


(0.20) 


o 


50 


(0.21) 


o 


41 


(0.23) 


o 


25 


(0.38) 


0.43 


68 


PMNJ0514-2029 


5 


321 


-5 


675 


2 


07 


(0.63) 


1 


43 


(0.82) 


1 


06 


(0.67) 


o 


41 


(0.45) 


o 


15 


(0.56) 


0.04 


71 


PMNJ0520-0537 


5 


427 


-48 


449 


1 


04 


(0.16) 


1 


25 


(0.24) 


1 


25 


(0.28) 


o 


63 


(0.25) 


1 


03 


(0.40) 


0.42 


74 


PMNJ0526-4830 


5 


680 


-54 


277 


1 


65 


(0.17) 


1 


52 


(0.20) 


1 


79 


(0.18) 


1 


22 


(0.21) 


o 


50 


(0.31) 


0.37 


77 


PMNJ0540-5418 


5 


839 


-57 


551 


1 


33 


(0.18) 


o 


94 


(0.21) 


1 


14 


(0.19) 


o 


54 


(0.21) 


o 


59 


(0.34) 


0.35 


79 


PMNJ0550-5732 


6 


003 


-45 


466 


o 


61 


(0.20) 


o 


95 


(0.21) 


o 


60 


(0.20) 


o 


66 


(0.23) 


o 


21 


(0.39) 


0.41 


81 


PMN T0559-4529 


8 


273 


-24 


423 


1 


11 


(0.18) 


o 


73 


(0.21) 


o 


87 


(0.21) 


o 


64 


(0.25) 


o 


32 


(0.40) 


0.27 


108 


PMN J0816- 2421 


10 


544 


41 


306 


1 


10 


(0.18) 


1 


02 


(0.20) 


o 


92 


(0.21) 


o 


87 


(0.26) 


o 


94 


(0.42) 


0.44 


133 


OB6 T1 033+41 15 

VJJJU O J- U '-J <J rt _1_ J_ ij 


11 


034 


-44 


Oil 


o 


70 


(0.19) 


o 


92 


(0.21) 


o 


77 


(0.20) 


o 


35 


(0.22) 


o 


23 


(0.34) 


0.45 


142 


PMNJ1102-4404 


11 


829 


-79 


554 


o 


97 


(0.18) 


o 


69 


(0.28) 


o 


51 


(0.24) 


o 


41 


(0.23) 


o 


10 


(0.41) 


0.04 


152 


PMN T1 1 50-791 8 


12 


056 


48 


104 


o 


79 


(0.24) 


o 


67 


(0.22) 


o 


67 


(0.23) 


o 


62 


(0.22) 


o 


14 


(0.34) 


0.16 


156 


GB6 T1 203+4803 


12 


464 


11 


406 





00 


(0.00) 





86 


(0.32) 





88 


(0.24) 





53 


(0.27) 





43 


(0.56) 


0.03 


162 


GB6 J1228+1124 


12 


519 


13 


857 





91 


(0.29) 





19 


(0.35) 





39 


(0.25) 





76 


(0.30) 





20 


(0.54) 


0.04 


165 


GB6 J1231+1344 


13 


047 


48 


947 





43 


(0.18) 





64 


(0.21) 





68 


(0.18) 





61 


(0.21) 


1 


09 


(0.31) 


0.03 


170 


GB6 J1303+4848 


13 


559 


27 


396 





88 


(0.21) 





81 


(0.26) 





68 


(0.27) 





52 


(0.25) 


1 


46 


(0.35) 


0.32 


181 


GB6 J1333+2725 


14 


668 


49 


973 





69 


(0.18) 


1 


01 


(0.22) 





84 


(0.20) 





45 


(0.25) 





74 


(0.33) 


0.20 


196 


GB6 J1439+4958 


16 


807 


41 


239 





80 


(0.22) 





71 


(0.31) 





71 


(0.27) 





29 


(0.23) 





00 


(0.38) 


0.20 


222 


GB6 J1648+4104 


16 


994 


68 


495 





42 


(0.17) 





64 


(0.18) 





72 


(0.16) 





79 


(0.19) 





52 


(0.29) 


0.38 


228 


GB6 J1700+6830 


17 


128 


1 


786 


1 


00 


(0.20) 





93 


(0.24) 





80 


(0.19) 





76 


(0.28) 





70 


(0.45) 


0.46 


230 


PMNJ1707+0148 


17 


618 


-79 


578 





82 


(0.18) 





78 


(0.20) 





81 


(0.18) 





74 


(0.21) 





69 


(0.38) 


0.42 


235 


PMNJ1733-7935 


22 


496 


-20 


853 





95 


(0.13) 





83 


(0.20) 





78 


(0.23) 





70 


(0.27) 


1 


28 


(0.47) 


0.33 


297 


PMNJ2229-2049 


22 


939 


-20 


191 


1 


00 


(0.20) 





61 


(0.22) 





76 


(0.24) 





80 


(0.30) 





41 


(0.45) 


0.45 


305 


PMNJ2256-2011 


23 


266 


-50 


286 





81 


(0.18) 


1 


12 


(0.18) 





77 


(0.19) 





80 


(0.23) 





61 


(0.35) 


0.34 


307 


PMNJ2315-5018 



NEWPS 6a SOURCE CATALOGUE 

The NEWPS 5(T Catalogue consists of 369 entries corresponding to all the EPS detected in the WMAP 3-yr full-sky 
maps at the > 5a level, after filtering with the MHW2 as discussed in the text. The 25 missed sources (WMAP detected 
sources that did not appear in our IC) are listed in Table 4 while the 369 IC sources detected at > 5<t are presented 
in Table 5. For each EPS in the Catalogue, and from left to right, we list the following data: the Equatorial, a, 5, 
coordinates of the center of the pixel; the NON CORRECTED source fluxes 7 , S v , in each WMAP channel identified 
by the channel symbols (K, Ka, Q, V and W) and their corresponding estimated errors (er), in Jy; the position in the 
3-year WMAP catalogue; the closest source present in the PMN or GB6 catalogues or the brightest NVSS or SUMSS 
source within the resolution element (i.e., inside the FWHM of the beam of the given WMAP frequency channel) and 
a "M" label that means that at least another source inside the beam has a 5 GHz flux within a factor of 2 from the 
brightest one, that we list as the likely identification. Table 4 has an additional column listing the 5 GHz fluxes. For 
those sources not detected at > 5ct, we have used "-" in the corresponding channel column. 



TABLE 5 

NEWPSscr SOURCE CATALOGUE 



RA 


Dec 




K(cr) 


Ka(<x) 


Q» 


V(<r) 


W(o-) 


WMAP 


low freq. id. 




h 


deg 


Jy(Jy) 


Jy(Jy) 


Jy(Jy) 


Jy(Jy) 


Jy(Jy) 








0.103 


-6.481 


2.7 


(0.2) 


2.3 (0.2) 


2.4 (0.2) 


2.2 (0.3) 




1 


PMNJ0006-0623 




0.214 


-39.957 


1.3 


(0.2) 


1.2 (0.2) 








2 


PMNJ0013-3954 




0.326 


26.031 


1.1 


(0.2) 










3 


NVSS J001939+260245 


M 


0.328 


20.445 


1.0 


(0.2) 










4 


GB6 J0019+2021 




0.424 


-26.067 


1.1 


(0.2) 










5 


PMNJ0025-2602 


M 


0.789 


-25.251 


1.2 


(0.2) 










9 


PMNJ0047-2517 




0.790 


-73.125 


1.9 


(0.2) 


1.4 (0.2) 


1.5 (0.2) 








PMNJ0047-7308 




0.827 


-57.637 


1.4 


(0.2) 


1.2 (0.2) 


0.9 (0.2) 






10 


PMNJ0050-5738 




0.847 


-9.429 


1.1 


(0.2) 










13 


PMNJ0050-0928 




0.854 


-6.828 


1.2 


(0.2) 










11 


PMNJ0051-0650 




0.962 


30.412 






1.2 (0.2) 










GB6 J0057+3021 




0.987 


-72.215 


2.3 


(0.2) 


1.4 (0.2) 


1.4 (0.2) 








PMNJ0059-7210 




1.099 


48.322 






1.2 (0.2) 










GB6 J0105+4819 




1.115 


-40.584 


1.8 


(0.1) 


1.9 (0.2) 


1.7 (0.2) 


1.1 (0.2) 




14 


PMNJ0106-4034 




1.140 


1.583 


2.6 


(0.2) 


2.3 (0.3) 


2.2 (0.2) 


2.0 (0.3) 




16 


PMNJ0108+0134 




1.152 


13.335 


1.4 


(0.2) 


1.1 (0.2) 








15 


GB6 J0108+1319 




1.270 


-11.686 


1.4 


(0.2) 


1.1 (0.2) 


1.3 (0.2) 


1.4 (0.3) 




18 


PMNJ0116-1136 




1.318 


-21.685 


0.8 


(0.2) 












PMNJ0118-2141 




1.361 


11.807 


1.2 


(0.2) 




1.1 (0.2) 






19 


GB6 J0121+1149 




1.379 


25.070 


1.0 


(0.2) 












GB6 J0122+2502 




1.421 


-0.082 


1.2 


(0.2) 










20 


PMNJ0125-0005 





7 The catalogue with the corrected fluxes can be obtained from max.ifca.unican.es/canicgo/NEWPS 
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TABLE 5 

NEWPS5a- SOURCE CATALOGUE 



1.544 


-16.863 


1.8 


(0.2) 


1.6 


CO 21 


1.8 


CO 21 


1.3 


CO 21 


1.551 


-52.009 


0.7 


(0.1) 














1.624 


47.842 


4.4 


(0.2) 


4.3 


(0.2) 


4.0 


CO 2) 


3.0 


CO 31 


1.626 


-24.473 


1.1 


(0.2) 


1.1 


(0.2) 


1.6 


CO 91 


1.4 


CO 31 


1.824 


5.973 


0.9 


(0.2) 














1.870 


22.043 


1.3 


(0.2) 


1.6 


(0.2) 










2.082 


-17.060 


0.8 


(0.2) 














2.083 


15.271 


1.6 


(0.2) 


1.5 


(0.2) 


1.3 


CO 2) 






2.083 


^9 99^ 


1.6 


(0.2) 


1.4 


(0.2) 


1.3 


CO 91 






2.183 


-51.009 


2.8 


(0.2) 


2.9 


(0.2) 


2.9 


CO 91 


2.6 


CO 21 


2.299 


73.874 


1.9 


(0.2) 


1.4 


(0.2) 


1.3 


CO 21 






2.317 


1.394 


1.0 


(0.2) 














2.355 


35.927 


1.2 


(0.2) 


1.1 


(0.2) 










2.380 


-34.726 


0.8 


(0.1) 














2.390 


43.035 


1.9 


(0.2) 


1.3 


(0.2) 


1.5 


CO 21 






2.529 


13.302 


1.4 


(0.2) 














2.627 


28 805 


3.5 


(0.2) 


2.9 


(0.2) 


3.3 


Co 31 


2.7 


CO 31 


2.644 


16.561 


1.5 


(0.2) 


1.5 


(0.3) 


1.5 


CO 21 


1.6 


CO 31 


2.666 


4.250 


0.9 


(0.2) 














2.858 


43.264 


1.1 


(0.2) 














2.890 


-54.680 


2.5 


(0.1) 


2.6 


(0.2) 


2.7 


CO 21 


2.1 


CO 21 


2.993 


-0.342 


1.1 


(0.2) 














3.065 


47.292 


1.3 


(0.2) 


1.4 


(0.2) 










3.068 


-62.246 


1.4 


(0.2) 


1.4 


(0.2) 


1.3 


CO 21 


1.2 


CO 21 


3.143 


4.099 


1.3 


(0.2) 














3.170 


-60.972 


1.1 


(0.2) 


1.3 


(0.2) 










3.196 


-76.889 


1.2 


(0.2) 


1.0 


(0.2) 










3.301 


41.831 










3.1 


CO 21 


2.0 


CO 3) 


3.329 


41.530 


11.4 


(0.2) 


8.3 


(0.2) 


7.1 


CO 21 


4.9 


CO 31 


3.361 


-37.023 


9.1 


(0.2) 


5.3 


(0.2) 


3.5 


CO 21 


1.7 


CO 21 


3.365 


12.347 


1.4 


(0.2) 














3.407 


-37.222 






4.8 


(0.2) 


2.1 


CO 21 






3.425 


22.397 


1.1 


(0.2) 














3.499 


-23 909 


1.0 


(0.2) 


1.2 


(0.2) 


1.1 


CO 21 






3.573 


-40.148 


1.6 


(0.2) 


1.4 


(0.2) 


1.4 


CO 21 


1.5 


CO 21 


3.604 


32 255 






1.9 


(0.3) 


2.1 


CO 31 


1.7 


CO 31 


3.606 


-13.048 


1.0 


(0.2) 






1.1 


CO 21 






3.654 


-1.767 


2.8 


(0.2) 


2.5 


(0.3) 


2.7 


CO 21 


1.5 


CO 31 


3.678 


-21.311 


1.2 


(0.1) 


1.1 


(0.2) 


1.1 


CO 21 


1.3 


CO 21 


3.813 


-27.772 


1.0 


(0.2) 










1.1 


CO 21 


3.985 


10.395 


1.4 


(0.2) 














4.069 


-36 080 


4.0 


(0.2) 


4.3 


(0.2) 


4.2 


CO 21 


3.9 


CO 21 


4.096 


-13.143 


1.9 


(0.2) 


2.0 


(0.2) 


1.6 


CO 21 


1.4 


CO 31 


4.111 


-38.449 






1.3 


(0.3) 










4.142 


-75.112 


0.8 


(0.2) 














4.173 


76.908 


1.1 


(0.2) 














4.302 


38.044 


3.8 


(0.4) 


3.4 


(0.4) 


2.7 


CO 31 


3.3 


CO 31 


4.392 


-1.300 


10.6 


(0.2) 


9.8 


(0.3) 


9.4 


CO 31 


8.2 


CO 31 


4.399 


41.839 


1.4 


(0.2) 


1.5 


(0.3) 










4.412 


-37.937 


1.5 


(0.2) 






1.5 


CO 21 


1.3 


CO 21 


4.416 


0.598 


1.9 


(0.2) 


1.9 


(0.2) 


1.8 


CO 31 






4.479 


-37.976 


1.2 


(0.2) 


1.5 


(0.2) 


1.2 


CO 21 


1.1 


CO 21 


4.555 


5.377 


2.6 


(0.2) 


2.8 


(0.2) 


2.6 


CO 21 


2.2 


CO 31 


4.615 


29.618 


3.5 


(0.2) 


2.5 


(0.2) 


1.8 


CO 21 


1.6 


CO 31 


4.633 


30 052 










1.3 


CO 21 






4.677 


-43.529 


3.0 


(0.2) 


2.5 


(0.3) 


2.4 


CO 21 


1.7 


CO 21 


4.709 


-0.250 


1.0 


(0.2) 






1.3 


CO 21 






4.820 


11.323 


2.2 


(0.2) 


2.3 


(0.2) 


2.3 


CO 21 


1.9 


CO 31 


4.832 


-80 985 


2.0 


(0.2) 


1.9 


(0.2) 


1.5 


CO 21 


1.6 


CO 21 


4.883 


-28.159 


1.7 


(0.2) 


1.6 


(0.2) 


1.3 


CO 21 


1.2 


CO 21 


4.925 


-46.276 


3.9 


(0.2) 


4.0 


(0.2) 


4.1 


CO 21 


3.8 


CO 21 


4.947 


-23.407 


2.8 


(0.2) 


2.7 


(0.2) 


2.7 


CO 21 


2.1 


CO 21 


5.016 


-2.023 


1.3 


(0.2) 


1.4 


(0.2) 


1.5 


CO 31 






5.114 


-61.220 


2.3 


(0.2) 


1.8 


(0.2) 


1.6 


CO 21 


1.2 


CO 21 


5.231 


-22.010 


0.9 


(0.2) 














5.257 


-45.940 


1.9 


(0.2) 














5.326 


-45.736 


6.4 


(0.3) 


5.0 


(0.2) 


4.0 


(0.2) 


2.8 


(0.2) 


5.384 


-36.484 


4.0 


(0.2) 


3.4 


(0.2) 


3.4 


(0.2) 


3.1 


(0.2) 


5.473 


21.500 


6.0 


(0.8) 


4.3 


(0.7) 










5.542 


7.487 










1.4 


(0.3) 






5.554 


48.382 


1.2 


(0.2) 


1.4 


(0.2) 










5.651 


-44.050 


5.6 


(0.2) 


5.5 


(0.2) 


5.4 


(0.2) 


4.3 


(0.2) 


5.710 


49.814 


1.7 


(0.2) 














5.929 


39.834 


2.8 


(0.2) 


1.9 


(0.2) 










6.124 


67.357 


0.9 


(0.2) 














6.126 


-6.362 


7.0 


(0.4) 


7.0 


(0.6) 


8.9 


(0.5) 


7.4 


(0.4) 


6.148 


-22.383 


0.9 


(0.2) 















2.6 (0.5) 



2.6 (0.3) 



4.0 (0.5) 



2.1 (0.4) 



2.2 (0.4) 



21 

22 
23 

24 

25 
26 
27 



28 
29 
30 
31 
32 
33 



34 



35 
36 
37 
38 

39 
40 



41 
42 
43 

44 
45 
46 
47 
48 
49 
50 
51 
52 
53 

54 

56 
57 
58 
59 



60 
61 

62 
63 
64 
65 
66 
67 
69 
70 
72 
73 



76 
78 
80 
82 

83 



PMNJ0132-1654 
PMNJ0133-5159 
GB6 J0136+4751 
PMNJ0137-2430 
PMNJ0149+0556 
GB6 J0152+2206 
PMNJ0204-1701 M 
GB6 J0204+1514 
GB6 J0205+3212 
PMNJ0210-5101 
GB6 J0217+7349 
PMNJ0219+0120 
GB6 J0221+3556 
PMNJ0222-3441 
GB6 J0223+4259 
GB6 J0231+1323 
GB6 J0237+2848 
GB6 J0238+1637 
PMNJ0239+0416 
GB6 J0251+4315 
PMNJ0253-5441 
PMNJ0259-0020 
GB6 J0303+4716 
PMNJ0303-6211 
PMNJ0308+0406 
PMNJ0309-6058 
PMNJ0311-7651 
GB6 J0318+4153 
GB6 J0319+4130 
PMNJ0321-3711 
GB6 J0321+1221 
PMNJ0324-3716 
GB6 J0325+2223 
PMNJ0329-2357 
PMNJ0334-4008 
GB6 J0336+3218 
PMNJ0336-1302 
PMNJ0339-0146 
PMNJ0340-2119 
PMNJ0348-2749 
GB6 J0358+1026 
PMNJ0403-3605 
PMNJ0405-1308 
PMNJ0406-3826 
PMNJ0408-7507 
NVSS J041045+765645 
GB6 J0418+3801 
PMNJ0423-0120 
GB6 J0423+4150 
PMNJ0424-3756 
PMNJ0423+0031 
PMNJ0428-3756 
PMNJ0433+0521 
GB6 J0437+2940 
GB6 J0438+3004 
PMNJ0440-4332 
PMNJ0442-0017 
GB6 J0449+1121 
PMNJ0450-8100 
PMNJ0453-2807 
PMNJ0455-4616 M 
PMNJ0457-2324 
PMNJ0501-0159 
PMNJ0506-6109 
PMNJ0513-2159 M 
PMNJ0515-4556 
PMNJ0519-4546a M 
PMNJ0522-3628 
NVSS J052830+213301 
PMNJ0532+0732 
GB6 J0533+4822 
PMNJ0538-4405 
GB6 J0542+4951 
GB6 J0555+3948 
GB6 J0607+6720 
PMNJ0607-0623 
PMNJ0608-2220 



M 



M 



M 



M 



M 



M 
M 



WMAP extragalactic point sources 



TABLE 5 

NEWPSsa SOURCE CATALOGUE 



6.162 


-15.727 


3.8 


(0.2) 


3.5 


CO 21 


3.4 


(0 2) 


2.3 


CO 2) 




84 


6.451 


-5.874 


1.4 


(0.2) 


















6.486 


-19.973 


1.4 


(0.2) 


1.2 


CO 21 


1.6 


(0 2) 


1.3 


CO 2) 




86 


6.585 


-75.258 


4.5 


(0.3) 


3.8 


(ft 3"! 


4.5 


(ft 2) 


3.2 


(ft 2"! 




88 


6.608 


-20.607 


1.2 


(0.2) 
















89 


6.653 


73.360 










1.1 


(ft 91 








90 


6.776 


44.907 


3.2 


(0.2) 


2.4 


(ft. A) 


2.1 


(ft 91 


1.7 


(ft 3) 
(ft. a) 




91 


6.839 


-16.564 


2.8 


(0.2) 


2.4 


(n 


2.1 


(ft 91 


1.6 


(ft 21 






7.350 


4.029 


1.0 


(0.2) 
















92 


7.369 


71.363 


1.8 


(0.2) 


1.8 


(ft 2) 


2.2 


(0 2) 


1.9 


CO 2) 




93 


7.418 


14.403 


1.0 


(0.2) 


















7.428 


-0.949 


1.2 


(0.2) 






1 3 


(ft 9) 








94 


7.566 


50 305 






1.3 


CO 31 












96 


7.635 


17.668 


1.4 


(0.2) 


1.4 


(ft 21 


1.3 


(ft 9"! 








97 


7.658 


1.627 


2.0 


(0.2) 


2.2 


CO 2) 


2.5 


(0 2) 


2.3 


Co 2) 




98 


7.684 


31.216 


1.3 


(0.2) 
















99 


7.719 


-67.474 


1.4 


(0.1) 
















100 


7.763 


10.217 


1.3 


(0.2) 
















101 


7.768 


-0.713 


1.3 


(0.2) 
















102 


7.802 


-36.071 


2.5 


(0.4) 


1.9 


CO 31 


1.6 


(0 2) 










7.847 


12.544 


2.7 


(0.2) 


2.3 


CO 3) 


2.6 


CO 31 


1.6 


CO 31 

\y.o) 




103 


7.886 


53.824 


1.1 


(0.2) 
















104 


7.953 


9.971 


1.7 


(0.2) 


1.7 


CO 2) 


1.6 


(0 3) 








105 


8.093 


-1.004 


1.0 


(0.2) 


















8.138 


-7.848 


1.4 


(0.2) 


1.4 


(ft 2") 


1.4 


(ft 21 


1.8 


(ft 31 
(u.o) 




106 


8.227 


48.188 


1.1 


(0.2) 
















107 


8.255 


36.557 


0.9 


(0.2) 


















8.411 


39 296 


1.3 


(0.2) 
















109 


8.428 


3.138 


1.4 


(0.2) 


1.5 


(ft ?1 


1.4 


(ft 3) 


1.4 


cn 3\ 




110 


8.512 


24.154 


1.7 


(0.2) 


1.4 


(0 3) 


1.6 


(0 3) 


2.0 


CO 31 




111 


8.614 


-20.320 


2.9 


(0.2) 


2.4 


(ft 2) 


2.1 


(ft 2) 


1.7 


(ft 31 




112 


8.682 


13.226 


1.9 


(0.2) 


2.1 


(ft 2) 


1.7 


(ft 3) 








114 


8.688 


70.908 


1.8 


(0.2) 


1.8 


(0 2) 


1.8 


(ft 2) 


1.6 


CO 21 




115 


8.917 


20.127 


3.9 


(0.2) 


4.3 


(0 2) 


3.9 


(ft 2) 


3.7 


CO 31 




116 


9.009 


-28.100 


1.2 


(0.2) 


















9.041 


-14.290 


1.2 


(0.2) 






1.2 


(ft 2) 
( V -^J 








117 


9.056 


46.835 


1.1 


(0.2) 


















9.078 


-57.536 


1.1 


(0.2) 


















9.132 


-20.302 


1.1 


(0.2) 
















118 


9.150 


1.320 


1.9 


(0.2) 


1.9 


(ft 2") 


1.8 


(ft 31 


1.5 


cn 31 




119 


9.157 


42.916 










1.3 


(ft 2) 








120 


9.299 


-12.131 


2.2 


(0.2) 


1.2 


(ft 2"! 


1.1 


(ft 2) 








122 


9.357 


44.628 


1.3 


(0.2) 


1.3 


(ft 2) 


1.2 


(ft 2) 








123 


9.362 


-26 360 


1.5 


(0.2) 


1.3 


(ft 2) 


1.1 


(ft 2) 








125 


9.365 


62.301 


0.9 


(0.1) 
















124 


9.378 


-39.992 


1.3 


(0.2) 


1.0 


(ft 2") 


1.0 


(ft 2) 










9.452 


39.010 


6.7 


(0.2) 


5.3 


(0 3) 


5.1 


(ft 2) 


4.0 


CO 31 




126 


9.816 


Aft fiKQ 


1.5 


(0.2) 


1.6 


(ft 91 

(ft. A) 


1.2 


(ft 91 








127 


9.939 


69.642 


1.4 


(0.2) 






1.0 


(ft 2) 








128 


9.954 


55.444 


1.1 


(0.2) 
















129 


9.976 


47.405 


1.7 


(0.2) 


1.4 


(ft 2) 


1.2 


(ft 2) 








130 


9.978 


-41.221 


0.8 


(0.2) 


















10.246 


-45.101 


1.0 


(0.2) 


















10.247 


23.040 


1.3 


(0.2) 
















131 


10.410 


-18.675 










1.1 


(0 2) 










10.620 


-29 623 


1.4 


(0.2) 


1.2 


(ft 2") 


1.1 


(ft 2) 








134 


10.645 


5.191 


2.0 


(0.2) 


2.0 


(ft 2) 


1.6 


(ft 2) 


1.5 


CO 31 




135 


10.692 


-47.705 


1.2 


(0.2) 
















137 


10.692 


6.185 


1.5 


(0.2) 


1.6 


(ft 2) 


1.5 


(ft 21 
( u -^) 








136 


10.798 


71.789 


1.1 


(0.2) 
















138 


10.804 


-19.122 


1.2 


(0.2) 
















139 


10.959 


-80.042 


2.0 


(0.2) 


2.1 


(ft 2) 


2.0 


(0 2) 
( U -^J 


2.1 


CO 21 




141 


10.971 


1.533 


4.6 


(0.2) 


4.3 


(0 2) 


4.6 


(ft 2) 


4.0 


CO 31 




140 


11.115 


-44.857 


1.6 


(0.2) 


1.5 


(ft 2) 


1.1 


(ft 2) 


1.3 


(ft 21 




143 


11.303 


-46.591 


0.9 


(0.2) 
















144 


11.314 


12.618 


1.1 


(0.2) 






1.1 


(ft 2) 








145 


11.448 


-18.972 


1.6 


(0.2) 


1.4 


(0.2) 


1.3 


(0.2) 








146 


11.500 


-14.735 


1.6 


(0.2) 


1.6 


(0.2) 


1.5 


(0.2) 








147 


11.514 


38.251 


1.5 


(0.2) 


1.1 


(0.2) 












148 


11.756 


-48.568 


0.7 


(0.1) 
















149 


11.780 


39.953 


0.9 


(0.2) 
















150 


11.785 


-38.167 


1.9 


(0.2) 


2.0 


(0.2) 


1.9 


(0.2) 


1.5 


(0.3) 




151 


11.890 


49.453 


2.0 


(0.2) 


1.7 


(0.2) 


1.9 


(0.2) 


1.6 


(0.2) 




153 


11.913 


81.082 


1.2 


(0.2) 
















154 


11.992 


29.230 


2.0 


(0.2) 


2.2 


(0.2) 


2.0 


(0.2) 


1.8 


(0.2) 




155 


12.149 


-24.126 


1.4 


(0.2) 
















157 


12.179 


-52.860 


2.6 


(0.2) 


1.6 


(0.2) 


1.1 


(0.2) 











PMNJ0609-1542 

PMNJ0627-0553 

PMNJ0629-1959 

PMNJ0635-7516 

PMNJ0636-2041 M 

GB6 J0639+7324 

GB6 J0646+4451 

PMNJ0650-1637 

PMNJ0721+0406 

GB6 J0721+7120 

GB6 J0725+1425 

PMNJ0725-0054 

GB6 J0733+5022 

GB6 J0738+1742 

PMNJ0739+0137 

GB6 J0741+3112 

PMNJ0743-6726 

GB6 J0745+1011 

PMNJ0745-0044 

PMNJ0748-3605 

GB6 J0750+1231 

GB6 J0753+5353 

PMNJ0757+0956 

NVSS J080537-005814 

PMNJ0808-0751 

GB6 J0813+4813 M 

GB6 J0815+3635 

GB6 J0824+3916 

PMNJ0825+0309 

GB6 J0830+2410 

PMNJ0836-2017 

GB6 J0840+1312 

GB6 J0841+7053 

GB6 J0854+2006 

PMNJ0900-2808 

PMNJ0902-1415 

GB6 J0903+4650 

PMNJ0904-5735 

PMNJ0906-2019 

PMNJ0909+0121 

GB6 J0909+4253 

PMNJ0918-1205 

GB6 J0920+4441 M 

PMNJ0921-2618 

GB6 J0921+6215 

PMNJ0922-3959 

GB6 J0927+3902 

GB6 J0948+4039 

GB6 J0955+6940 

GB6 J0957+5522 

GB6 J0958+4725 

PMNJ0958-4110 

PMNJ1014-4508 

GB6 J1014+2301 

PMNJ1024-1838 

PMNJ1037-2934 

PMNJ1038+0512 

PMNJ1041-4740 

PMNJ1041+0610 

GB6 J1048+7143 

PMNJ1048-1909 

PMNJ1058-8003 

PMNJ1058+0133 

PMNJ1107-4449 M 

PMNJ1118-4634 

GB6 J1118+1234 

PMNJ1127-1857 

PMNJ1130-1449 

GB6 J1130+3815 

PMNJ1145-4836 

GB6 J1146+3958 

PMNJ1147-3812 M 

GB6 J1153+4931 

NVSS J115312+805829 M 

GB6 J1159+2914 

PMNJ1209-2406 

PMNJ1211-5250 



Lopez-Caniego et al. 



TABLE 5 

NEWPSsa SOURCE CATALOGUE 



5.5 (0.5) 
4.2 (0.4) 



12.262 


-17.437 


1.7 


(0.2) 


1.4 


CO 21 










12.322 


5.733 


2.4 


(0.2) 


2.2 


CO 3) 


2.0 (0.3) 


1 


5 


CO 31 


12.484 


2.085 


18.7 


(0.3) 


16.4 


CO 31 


15 2 (ft 31 


12 


3 


(ft 31 


12.507 


12.055 


19.3 


(0.3) 


14.7 


(ft 3"! 


1 2 2 (ft 21 


8 


1 


(ft 31 


12.762 


-16.280 


0.9 


(0.2) 














12.781 


-25.780 


1.4 


(0.2) 


1.4 


CO 21 


1.8 (0.2) 








12.937 


-5.742 


16.7 


(0.3) 


16.1 


(ft 3"! 


17.4 (0.3) 


14 


7 


CO 31 


12.966 


-31.981 


1.5 


(0.1) 


1.1 


CO 2) 


1.3 (0.2) 








12.979 


-22 320 


1.0 


(0.2) 














13.159 


11.942 


1.0 


(0.2) 














13.176 


32.356 


3.2 


(0.2) 


3.1 


(0 2) 


2.8 (0.2) 


2 


4 


CO 21 


13.266 


-33.624 


1.6 


(0.2) 


1.5 


(ft 2) 


1.6 (0.2) 


1 


3 


(ft 21 


13.351 


-43.707 






2.4 


(ft d\ 










13.388 


-44.845 


3.5 


(0.4) 


1.9 


(ft 41 










13.428 


-43.042 


46.2 


(0.5) 


36.7 


(ft 41 


34 8 (ft 31 


25 


2 


CO 31 


13.447 


-52.981 


1.8 


(0.2) 


1.8 


(ft 2) 


1.6 (0.2) 


1 


8 


CO 21 


13.453 


-42.678 










7.3 (0.3) 








13.457 


22 202 


1.1 


(0.2) 














13.499 


31.859 


1.0 


(0.2) 














13.516 


25.071 


1.1 


(0.2) 














13.518 


30.570 


2.3 


(0.2) 


1.9 


(ft 2) 


1.4 (0.2) 








13.552 


2.020 


1.3 


(0.2) 














13.624 


-34.130 


1.9 


(0.2) 


1.2 


(0.2) 










13.628 


-12.897 


6.2 


(0.2) 


6.2 


(0.2) 


6.1 (0.2) 


5 


3 


CO 31 


13.743 


66 066 

UU- UUU 


0.9 


(0.2) 














13.794 


12.303 


0.9 


(0.2) 


1.1 


(0.2) 










13.912 


-10.653 


1.7 


(0.2) 


1.3 


(0.2) 


1.7 (0.2) 








13.951 


76.746 


0.8 


(0.2) 














13.952 


19.315 


1.4 


(0.2) 


1.4 


(0.2) 


1.1 (0.2) 








13.952 


-15.488 


1.1 


(0.2) 














14.148 


-7.826 


1.3 


(0.2) 


1.1 


(0.2) 


1.2 (0.2) 








14.268 


13.325 


0.9 


(0.2) 














14.327 


38 385 


0.9 


(0.2) 


1.0 


(0.2) 




1 





(ft 21 


14.327 


54.355 


1.0 


(0.2) 








1 


2 


(ft 21 


14.334 


27.040 


1.1 


(0.1) 














14.408 


-49.221 


2.2 


(0.2) 


1.6 


(0.2) 


1.5 (0.2) 








14.412 


-68.110 


1.1 


(0.2) 






1.2 (0.2) 








14.457 


-33.351 


1.0 


(0.2) 


1.5 


(0.2) 


1.8 (0.2) 


1 


5 


CO 31 


14.463 


-42.116 


3.2 


(0.2) 


2.9 


(0.2) 


2 8 (ft 2*1 


2 


3 


CO 31 


14.635 


-22.120 


0.8 


(0.2) 


1.1 


(0.2) 










14.714 


52 080 


1.0 


(0.2) 














14.909 


-37.741 


1.0 


(0.2) 














14.991 


71.762 


1.5 


(0.2) 


1.6 


(0.2) 


1.0 (0.2) 








15.064 


-41.927 


2.3 


(0.2) 


1.6 


(0.2) 


1.6 (0.2) 








15.078 


10.490 


2.0 


(0.2) 


1.7 


(0.2) 


1.5 (0.2) 








15.119 


-16.887 


1.6 


(0.2) 


1.5 


(0.2) 










15.180 


-5.634 


1.1 


(0.2) 














15.211 


-9.147 


1.8 


(0.2) 


1.7 


(0.2) 


1.8 (0.2) 


1 


7 


CO 3"! 


15.228 


-10.228 


1.0 


(0.2) 














15.248 


-47.857 


1.7 


(0.2) 


1.4 


(0.2) 


1.4 (0.2) 








15.285 


0.206 


1.7 


(0.2) 


1.8 


(0.2) 


1.6 (0.2) 








15.296 


-24.422 


2.0 


(0.2) 


1.9 


(0.2) 


1.9 (0.2) 


1 


p 


(ft 31 


15.380 


-27.510 


1.3 


(0.3) 






1.2 (0.2) 








15.680 


14.846 


1.2 


(0.2) 














15.825 


50.657 


1.1 


(0.2) 






9 (ft 2) 








15.826 


2.660 


2.3 


(0.2) 


2.3 


(0.3) 


2 (ft 21 


1 


g 


(ft 31 


15.845 


5.469 


2.6 


(0.2) 


2.2 


(0.3) 


1.8 (0.3) 


1 


9 


CO 31 


16.042 


33.446 


0.9 


(0.1) 














16.150 


10.456 


2.6 


(0.2) 


2.4 


(0.2) 


2 (ft 2*1 


1 


7 


CO 31 


16.228 


34.206 


3.9 


(0.2) 


3.3 


(0.2) 


3 2 (ft 2) 


2 


7 


CO 21 


16.250 


-60.657 


2.6 


(0.2) 


1.8 


(0.3) 










16.286 


-58.815 


2.6 


(0.4) 






2.1 (0.4) 








16.310 


-77.341 


2.4 


(0.2) 


2.2 


(0.2) 


1.9 (0.2) 


1 


7 


CO 21 


16.431 


-25.425 










2 3 (ft 31 


2 


1 


CO 31 


16.431 


-29.903 


1.6 


(0.3) 














16.539 


82 520 


1.2 


(0.1) 


1.3 


(0.2) 


1.2 (0.2) 








16.581 


38.150 


4.2 


(0.2) 


4.7 


(0.2) 


4.5 (0.2) 


3 


6 


(0.2) 


16.636 


47.260 










1.0 (0.2) 








16.640 


57.362 


1.3 


(0.2) 


1.4 


(0.2) 


1.5 (0.2) 


1 


4 


(0.2) 


16.676 


39.774 










1.8 (0.2) 








16.708 


68.891 


1.2 


(0.2) 


1.0 


(0.2) 


1.2 (0.2) 


1 


3 


(0.2) 


16.719 


39.812 


7.0 


(0.3) 


6.2 


(0.2) 


6.0 (0.2) 


5 


2 


(0.3) 


16.855 


4.951 


1.9 


(0.2) 


1.3 


(0.2) 


1.2 (0.2) 








16.897 


39.715 


1.2 


(0.2) 














16.965 


47.823 


0.9 


(0.2) 














16.969 


7.628 


1.1 


(0.2) 














17.026 


-56.367 


1.2 


(0.2) 















7.9 (0.5) 



11.5 (0.4) 



2.2 (0.4) 



2.7 (0.4) 



158 
160 
163 
164 

166 
167 
168 
169 
172 
173 
174 



176 
177 
178 
179 
180 
182 
183 
184 

185 
187 
186 
188 
189 

192 
191 
193 



194 
195 

197 

198 

199 
200 
201 
202 
203 

204 
205 

206 
207 
208 
209 
211 
212 
213 



214 



215 
216 
217 
218 

219 
220 
223 
224 
226 
227 



PMNJ1215-1731 
PMNJ1219+0549 
PMNJ1229+0203 
GB6 J1230+1223 M 
PMNJ1245-1616 
PMNJ1246-2547 
PMNJ1256-0547 
PMNJ1257-3154 
PMNJ1258-2219 
GB6 J1309+1154 
GB6 J1310+3220 M 
PMNJ1316-3339 
PMNJ1321-4342 
PMNJ1323-4452 
PMNJ1325-4257 M 
PMNJ1326-5256 
PMNJ1327-4239 
GB6 J1327+2210 
GB6 J1329+3154 
GB6 J1330+2509 
GB6 J1331+3030 
GB6 J1332+0200 
PMNJ1336-3358 M 
PMNJ1337-1257 
GB6 J1344+6606 
GB6 J1347+1217 
PMNJ1354-1041 
NVSS J135755+764320 
GB6 J1357+1919 
PMNJ1357-1527 
NVSS J140856-075226 
GB6 J1415+1320 
GB6 J1419+3822 M 
GB6 J1419+5423 
NVSS J141958+270143 
PMNJ1424-4913 
PMNJ1424-6808 
NVSS J142717-331800 
PMNJ1427-4206 
PMNJ1438-2204 
GB6 J1443+5201 
PMNJ1454-3747 
GB6 J1459+7140 
PMNJ1503-4154 
GB6 J1504+1029 
PMNJ1507-1652 
PMNJ1510-0543 
NVSS J151250-090600 
PMNJ1513-1012 
PMNJ1514-4748 
PMNJ1516+0014 
PMNJ1517-2422 
PMNJ1522-2730 
GB6 J1540+1447 
GB6 J1549+5038 
PMNJ1549+0237 
PMNJ1550+0527 
GB6 J1602+3326 
GB6 J1608+1029 
GB6 J1613+3412 
PMNJ1615-6054 
PMNJ1617-5848 
PMNJ1617-7717 
PMNJ1625-2527 
NVSS J162605-295126 
NVSS J163051+823345 
GB6 J1635+3808 
GB6 J1637+4717 
GB6 J1638+5720 
GB6 J1640+3947 
GB6 J1642+6856 
GB6 J1642+3948 
PMNJ1651+0459 
GB6 J1653+3945 M 
GB6 J1658+4737 
PMNJ1658+0741 
PMNJ1701-5621 



M 
M 



M 



WMAP extragalactic point sources 



TABLE 5 

NEWPSsa SOURCE CATALOGUE 



17.064 


-62 252 


1.8 


(0.2) 


1.8 


CO 2"! 


1.9 


(ft 21 


1.4 


CO 21 


17.341 


-1.005 


5.7 


(0.2) 


3.9 


(0 2) 


3.3 


(0 2) 


2.6 


(0 31 


17.388 


-65 003 


2.4 


(0.2) 


1.9 


CO 21 


1.5 


(ft 21 


1.2 


CO 21 


17.453 


45.540 


1.1 


(0.2) 


1 3 


(ft 21 






1.3 


(ft 21 


17.488 


-23.821 










2.1 


(ft 41 






17.515 


-21.458 


2.2 


(0.4) 


1.7 


(0 3) 


1.9 


(0 4) 






17.551 


-13.075 


5.4 


(0.2) 


5.4 


(ft 21 


5.1 


CO 21 


4.0 


en 31 


17.570 


39.009 


1.1 


(0.2) 


1.3 


CO 21 


1.1 


CO 21 


1.2 


(ft 21 


17.625 


6.316 


1.3 


(0.2) 






1.1 


(ft 21 






17.663 


-71.251 










0.9 


(0 21 






17.667 


47.672 










1.0 


(ft 21 






17.678 


52.175 


1.3 


(0.2) 


1.1 


CO 2) 


1.3 


(ft 21 


1.4 


CO 21 


17.731 


-3.789 


5.8 


(0.2) 


5.2 


(ft 2) 


4.9 


CO 21 


4.5 


CO 31 


17.739 


-51.762 


1.0 


(0.2) 














17.797 


70.113 














1.0 


CO 21 


17.860 


9.709 


4.1 


(0.2) 


4.2 


(0 2) 


4.3 


(0 21 


4.2 


CO 21 


17.879 


9.941 














2.0 


(ft 21 


17.884 


9.351 










1.4 


(ft 21 






17.888 


44.131 


0.9 


(0.2) 














17.893 


28.730 


2.3 


(0.1) 


2.2 


(ft 21 


2.3 


CO 21 


2.3 


CO 21 


18.008 


38.842 


1.0 


(0.1) 














18.014 


78.470 


2.1 


(0.1) 


1.6 


(ft 2) 


1.7 


(0 21 


1.5 


CO 21 


18.023 


44.067 


1.1 


(0.2) 


1.1 


(0 2) 


1.4 


(0 2) 


1.4 


CO 21 


18.047 


-39.695 


1.9 


(0.2) 


1.9 


(ft 2} 


1.6 


(ft 21 


2.0 


(ft 31 


18.056 


-65.156 


1.3 


(0.2) 


1.2 


(ft 2) 


1.5 


(ft 21 






18.119 


69.806 


1.6 


(0.1) 


1.5 


(0 2) 


1.3 


(0 2) 


1.5 


CO 21 


18.166 


-45.918 


1.5 


(0.2) 


1.5 


(ft 21 


1.3 


(ft 21 


1.3 


CO 31 


18.205 


6.838 


1.1 


(0.2) 














18.324 


-55.350 


0.9 


(0.2) 














18.363 


-63 968 


1.6 


(0.2) 


1.5 


(ft 21 










18.398 


56.812 


1.5 


(0.2) 


1.4 


(ft 21 


1.4 


CO 21 


1.2 


CO 21 


18.495 


48.771 


2.7 


(0.2) 


2.4 


(0 2) 


2.4 


(0 2) 


1.6 


CO 21 


18.562 


-21.110 


4.9 


(0.5) 


4.5 


(0 51 


4.6 


CO 51 


2.7 


CO 41 


18.584 


32 650 


0.9 


(0.1) 














18.609 


-71.853 


0.8 


(0.1) 














18.626 


-71.134 


1.8 


(0.1) 


1.5 


(0 2) 


1.3 


(0 21 


1.0 


CO 21 


18.719 


79.780 


1.1 


(0.2) 














18.722 


68.146 


1.0 


(0.2) 


1.2 


(ft 2) 


1.1 


(0 21 






18.824 


67.155 


1.4 


(0.2) 


1.5 


(ft 2) 


1.3 


("0 21 


1.0 


CO 21 


18.842 


28 396 


1.2 


(0.2) 














19.051 


31.972 


1.0 


(0.2) 














19.183 


-20.091 


2.0 


(0.2) 


1.9 


(ft 31 


1.7 


(0 31 


1.7 


CO 31 


19.394 


-21.033 


2.3 


(0.3) 


2.1 


(0 31 


2.0 


(ft 31 


2.0 


CO 31 
yv-o) 


19.417 


-29.201 


12.1 


(0.2) 


10.9 


(ft 31 


10.6 


en 21 


9.6 


(ft 31 


19.455 


61.295 


1.2 


(0.2) 


1.1 


(0 21 


1.2 


Co 21 






19.460 


73.985 


3.6 


(0.2) 


3.2 


(0 21 


2.6 


(0 2) 


2.5 


CO 21 


19.621 


-40.005 










1.2 


(ft 2) 






19.657 


-15.422 


1.1 


(0.2) 














19.659 


-63.710 


0.9 


(0.2) 














19.659 


-69.144 


0.7 


(0.1) 














19.872 


2.514 


1.1 


(0.2) 














19.934 


51.548 






1.1 


(ft 21 
1 u - z 1 










19.961 


-38.736 


3.3 


(0.2) 


3.2 


(ft 21 


2.9 


(ft 21 


2.5 


CO 31 


19.964 


-55.167 


0.9 


(0.2) 














19.997 


40.704 


54.6 


(0.4) 


34.9 


(0 41 


30.4 


("0 41 


18.0 


CO 41 


20.018 


-17.850 


1.8 


(0.2) 


1.5 


(ft 21 


1.3 


CO 21 


1.5 


CO 31 


20.100 


77.924 






1.0 


(ft 21 


0.9 


co 21 






20.159 


-48.794 


0.8 


(0.2) 














20.190 


-15.743 


1.6 


(0.2) 


1.5 


(0 31 


1.5 


(0 21 






20.372 


61.679 


1.6 


(0.2) 


1.4 


(ft 21 


1.1 


(ft 21 






20.415 


17.256 


0.9 


(0.2) 














20.609 


-68.771 


0.9 


(0.1) 


1.1 


(0 21 










20.640 


51.296 


2.4 


(0.3) 






1.5 


(0 21 


1.7 


CO 31 


20 933 


-47.226 


1.6 


(0.2) 


1.5 


CO 21 
1 u - z 1 


1.7 


(ft 21 


1.5 


CO 31 


21.024 


3.682 


0.8 


(0.2) 














21.034 


-27.992 


1.0 


(0.2) 














21.159 


35.494 


1.2 


(0.2) 














21.166 


-41.208 


1.3 


(0.2) 


1.6 


(0-2) 


1.2 


(0.2) 






21.394 


5.611 


2.2 


(0.2) 


1.8 


(0.2) 


1.7 


(0.2) 


1.3 


(0.3) 


21.396 


25.014 


0.9 


(0.2) 














21.529 


-12.153 


2.9 


(0.2) 


2.3 


(0-2) 


2.1 


(0.2) 


1.8 


(0.3) 


21.570 


-1.904 


1.9 


(0.2) 


1.8 


(0-2) 


1.7 


(0.3) 


1.5 


(0.3) 


21.606 


0.734 


4.5 


(0.3) 


3.7 


(0-3) 


2.9 


(0.3) 






21.651 


14.410 


2.2 


(0.2) 


1.9 


(0-2) 


1.6 


(0.2) 






21.730 


17.741 


1.2 


(0.2) 


1.1 


(0.2) 










21.762 


-77.935 


1.3 


(0.2) 


1.1 


(0.2) 










21.803 


7.015 


8.1 


(0.2) 


7.2 


(0.2) 


6.9 


(0.2) 


5.6 


(0.3) 



2.5 (0.4) 



4.4 (0.4) 



7.9 (0.5) 



3.0 (0.4) 



229 



232 
233 



234 



236 
237 



238 



239 
240 
242 
243 
244 

245 
246 



247 
248 
249 
250 

251 

252 
253 
254 
256 
257 
258 

259 

260 
261 
262 

263 

264 
265 
266 



267 
268 

270 
271 

272 

274 



276 
275 
277 

278 
279 
280 
281 
282 
284 
283 



PMNJ1703-6212 M 

PMNJ1720-0058 M 

PMNJ1723-6500 

GB6 J1727+4530 M 

PMNJ1729-2345 

PMNJ1730-2129 

PMNJ1733-1304 

GB6 J1734+3857 

GB6 J1737+0620 

SUMSS J173939-712325 

GB6 J1739+4738 

GB6 J1740+5211 5 

PMNJ1743-0350 

PMNJ1744-5144 

GB6 J1748+7005 

GB6 J1751+0938 M 

NVSS J175231+095711 

NVSS J175302+091959 

GB6 J1753+4410 M 

GB6 J1753+2847 

GB6 J1800+3848 M 

NVSS J180045+782805 

GB6 J1801+4404 

PMNJ1802-3940 

PMNJ1803-6507 M 

GB6 J1806+6949 

PMNJ1809-4552 

GB6 J1812+0651 

PMNJ1819-5521 

PMNJ1819-6345 

GB6 J1824+5650 

GB6 J1829+4844 

PMNJ1833-2103 

GB6 J1835+3241 

PMNJ1835-7150 

PMNJ1837-7108 

NVSS J184226+794517 

GB6 J1842+6809 

GB6 J1849+6705 

GB6 J1850+2825 

GB6 J1902+3159 

PMNJ1911-2006 

PMNJ1923-2104 

PMNJ1924-2914 M 

GB6 J1927+6117 

GB6 J1927+7357 

PMNJ1937-3957 

PMNJ1939-1525 M 

PMNJ1939-6342 

PMNJ1940-6908 

PMNJ1952+0230 

GB6 J1955+5131 

PMNJ1957-3845 

PMNJ1958-5509 

NVSS J195921+403428 M 

PMNJ2000-1748 

NVSS J200531+775243 

PMNJ2009-4849 

PMNJ2011-1546 

GB6 J2022+6137 

GB6 J2024+1718 

PMNJ2035-6846 

GB6 J2038+5119 

PMNJ2056-4714 

PMNJ2101+0341 

PMNJ2101-2802 

GB6 J2109+3532 

PMNJ2109-4110 

PMNJ2123+0535 

GB6 J2123+2504 

PMNJ2131-1207 

PMNJ2134-0153 

PMNJ2136+0041 

GB6 J2139+1423 

GB6 J2143+1743 M 

PMNJ2146-7755 

PMNJ2148+0657 
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NEWPSsa SOURCE CATALOGUE 



21.855 


7.113 


1.3 


(0.2) 


















PMN T21 51 +0709 


21.872 


-30.498 


1.1 


(0.2) 


1.3 


en 2i 


1.2 


CO 21 


1.5 


CO 31 




285 


PMN 121 51-3028 

x ±vx l^ jiiiui oUiiu 


21.888 


47.306 


2.5 


(0.2) 


2.7 


en 3i 


2.2 


CO 21 


2.1 


CO 21 


2.0 (0.4) 




OR6 12153+4716 


21.962 


-69 639 


3.8 


(0.2) 


2.8 


en 21 


2.6 


CO 21 


1.9 


CO 21 




286 


PMN 121 57-6941 


21.969 


-15.057 


2.3 


(0.2) 


2.0 


en 21 


1.8 


CO 21 








287 


PMN 121 58-1 501 

X 1VX1 N O irf_LOO 1UU1 


22.048 


42.317 


3.1 


(0.2) 


2.8 


CO 21 


2.9 


CO 21 


2.4 


CO 21 




288 


GB6 J2202+4216 


22 053 


31.777 


2.9 


(0.2) 


2.8 


en 21 


2.1 


CO 21 


1.7 


CO 21 




289 


CB6 12203+3145 


22.059 


17.410 


1.7 


(0.2) 


1.6 


en 21 


1.6 


CO 21 


1.7 


CO 31 




290 


GB6 12203+1725 


22.101 


-18.642 


1.8 


(0.2) 


1.5 


en 21 


1.2 


CO 21 








291 


PMN 12206-1 835 


22.133 


-53.712 


1.1 


(0.2) 
















292 




22.199 


23.874 


1.1 


(0.2) 


1.3 


CO 21 












293 


C1B6 12212+2355 


22.218 


-25.529 


0.9 


(0.2) 


















PMN 1221 3-2529 


22.314 


-3.630 


2.7 


(0.2) 


2.0 


CO 21 


2.0 


CO 41 








294 


PMN 1221 8-0335 


22.426 


21.289 






1.2 


CO 21 














GB6 J2225+2118 


22.427 


-4.985 


5.0 


(0.2) 


4.6 


en 31 


3.9 


CO 31 


3.5 


CO 31 




295 


PMN 12225-0457 


22.498 


-8.534 


2.5 


(0.2) 


2.6 


en 31 


2.7 


CO 31 


3.0 


CO 31 




296 


PMN 12229-0832 


22.547 


11.767 


3.4 


(0.2) 


3.2 


CO 21 


3.4 


CO 21 


3.1 


CO 31 


2 6 CO 51 


298 


GB6 J2232+1143 


22.589 


-48.650 


1.9 


(0.2) 


1.8 


CO 21 


1.7 


CO 21 


1.5 


CO 21 




299 


PMNJ2235-4835 


22.603 


28.502 


1.3 


(0.2) 


1.3 


CO 21 












300 


GB6 J2236+2828 


22.648 


-56.994 


0.9 


(0.2) 
















301 


PMN 12239-5701 


22.774 


-12.133 


1.2 


(0.2) 


1.2 


CO 21 


1.1 


CO 21 








302 


PMN J2246- 1206 


22.901 


16.135 


7.4 


(0.2) 


6.7 


CO 21 


6.8 


CO 21 


6.1 


CO 31 


4.0 (0.4) 


303 


OB6 12253+1608 


22.958 


41.923 


1.0 


(0.2) 


















GB6 J2257+4154 


22.970 


-28 003 


6.8 


(0.2) 


6.4 


CO 21 


6.2 


CO 21 


5.6 


CO 31 


3.5 (0.4) 


306 


PMN 12258-2758 


23.367 


27.500 


0.9 


(0.2) 


















GB6 J2322+2732 


23.395 


-3.296 


1.0 


(0.2) 


















PMN 12323-031 7 


23.457 


9.663 






1.4 


CO 21 


1.3 


CO 21 








309 


PMNJ2327+0940 


23.489 


-47.438 


1.3 


(0.2) 


1.1 


CO 21 


1.3 


CO 21 








310 


PMN 12329-4730 


23.512 


10.973 


1.0 


(0.2) 
















311 


vjDU J ZOOUT-L 1UU 


23.528 


-15.973 


1.4 


(0.2) 




_ 




_ 




_ 


_ 


312 


PMNJ2331-1556 


23.566 


-23.721 










1.2 


(0.2) 








313 


PMNJ2333-2343 


23.568 


7.559 


1.2 


(0.2) 
















314 


PMNJ2334+0736 


23.589 


-1.493 






1.2 


(0.2) 


1.3 


(0.2) 








315 


PMNJ2335-0131 


23.602 


-52.602 


0.9 


(0.2) 
















316 


PMNJ2336-5236 


23.778 


9.453 


1.1 


(0.2) 
















317 


PMNJ2346+0930 


23.803 


-16.538 


1.9 


(0.2) 


1.7 


(0.3) 


2.0 


(0.2) 


1.6 


(0.3) 




318 


PMNJ2348-1631 


23.907 


45.943 


1.7 


(0.1) 


1.2 


(0.2) 


1.2 


(0.2) 








319 


GB6 J2354+4553 


23.964 


-53.162 


1.3 


(0.1) 


1.1 


(0.2) 


1.2 


(0.2) 








321 


PMNJ2357-5311 


23.964 


-10.295 


1.0 


(0.2) 
















322 


PMNJ2358-1020 


23.979 


-60.967 


1.8 


(0.2) 


1.4 


(0.2) 


1.1 


(0.2) 








323 


PMNJ2358-6054 



